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COREST: A FORTRAN COMPUTER PROGRAM TO ANALYZE 
PARALINEAR OXIDATION BEHAVIOR AND ITS APPLICATION 
TO CHROMIC OXIDE FORMING ALLOYS 
by Charles A. Barrett and Alden F. Presler 
Lewis Research Center 

SUMMARY 

A FORTRAN program termed "COREST" has been developed for analyzing para- 
linear oxidation. It is based on the mass-balance expression AW (net sample weight 
change) = W r - W m . The W r (specific weight of the retained oxide) and W m (specific 
weight of the metal consumed) values can be generated independently as a function of 
time from a minimal amount of input weight change data, assuming that the type of oxide 
which forms is known. The W m value is critical in corrosion. It is always increas- 
ingly greater with time in paralinear oxidation than in parabolic oxidation with the same 
parabolic scaling constant. 

This technique of paralinear analysis was evaluated and demonstrated for predomi- 
nantly CrgOg-forming alloys tested for long times above 1100° C. In general, experi- 
mental AW- against- time data agree quite closely with values estimated by the COREST 
paralinear program, usually within ±0. 5 mg/cm . The derived paralinear scaling and 
vaporization rate constants kp and k y values were then used to predict long-time 
metal consumption for Cr 202 /chromite- forming alloys in the 1000° to 1200° C range. 
Comparisons with parabolic oxidation data analysis show that the assumption of para- 
linear kinetics is more realistic for the alloys examined, particularly for long-time 
predictions. The paralinear model includes the parabolic, where the vaporization con- 
stant k y approaches zero, as the limiting case. It is suggested that paralinear analysis 
be performed first on shorter time data to indicate how long a test should be run to 
meaningfully describe the oxidation behavior. 


INTRODUCTION 


The high- temperature, isothermal oxidation of metals and alloys is generally ana- 



lyzed on the basis of oxygen weight gain as a function of time. This process follows par- 
abolic kinetics if the oxide scale grows under conditions which are controlled by bulk dif- 
fusion. In this case, metal consumption with time can be calculated directly. In many 
cases, however, parabolic kinetics are not followed. For many high- chromium (Cr)- 
content, heat-resistant iron (Fe), nickel (Ni), or cobalt (Co) alloys, the scale formed is 
predominantly CrgOg. This oxide can further react with oxygen and volatilize according 
to the reaction CrgO^s) + Og(g) — 2CrOg(g) (ref. 1). This volatilization process is of 
special concern in the high-velocity environments that are encountered by orbital reen- 
try vehicles and gas turbine engines. Thus, weight gains due to oxygen pickup and 
weight losses due to CrOg volatilization occur simultaneously. In this situation, the 
sum of the weight gain and weight loss processes does not produce a weight change curve 
that is parabolic with time. Therefore, metal loss (i.e. , corrosion) cannot be directly 
estimated. Here, typical weight- change -against- time curves first increase to a maxi- 
mum, then show a gradual weight loss, and finally exhibit an almost linear weight loss. 
Such behavior belongs to a general class of reactions that are termed "paralinear. M 

The paralinear reaction was defined by Haycock (ref. 2) for a barrier scale layer 
being formed by diffusion and simultaneously being depleted at a constant rate by a sec- 
ondary process. Haycock derived a general solution for the rate equation, and Wajszel 
(ref. 3) solved it specifically for a diffusion- controlled, volatile -scale -formation case. 
More recently, beginning with Tedmon (ref. 4), Giggins and Pettit (ref. 5), Lewis 
(ref. 6), Stringer (ref. 7), and Ecer and Meiser (ref. 8), the problem of analyzing para- 
linear data and deriving the proper constants has been investigated. Because of the 
nature of the equations, the experimental data are difficult to interpret and analyze un- 
less a convenient computer program is available. 

This report presents a computer program that can be used to analyze and extrapo- 
late paralinear oxidation data in a convenient, readily usable form. This program, 
called COREST, is flexible enough so that it can be used under various conditions of ox- 
idation testing and analysis and still estimate corrosion attack. In this report, para- 
linear oxidation analysis and the COREST computer program are described in detail. 

The COREST program is then applied to the analysis of the long-time oxidation test re- 
sults for various CrgOg/chromite-forming alloys that were available in the literature or 
that were obtained at this laboratory. 

Included is a comparison of how well the calculated specific weight change values 
match the experimental data points depending on whether short-, middle-, or long-time 
input data are used. The comparison is intended to answer the critical questions - how 
long should an oxidation test be run to properly interpret oxidation behavior and how can 
one make the distinction between parabolic and paralinear kinetics ? 

Finally, the derived rate constants are used in the COREST program to estimate the 
long-time (100 000 hr) oxidation behavior of CrgOg/ chromite-forming alloys in the 1000° 
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to 1200° C range, with emphasis on metal consumption. These estimates assume that 
paralinear kinetics still control oxidation at times beyond 10 000 hours. These esti- 
mates were made for the most and the least oxidation resistant of the predominantly 
Cr 2®3 _forming all °y g considered for several airflow (i.e. , vaporization) conditions. 
Long-time oxidation behavior estimates were also made by assuming parabolic kinetics, 
and the results were compared with those obtained from the COREST program. 


SYMBOLS 


A 


a 


B 

b 


parabolic scaling constant for oxidation in terms of weight of metal consumed 

n o 

(not the usual weight of oxygen picked up), (mg/cm ) /hr (So in terms of the 
rate constant kp in the oxidation literature, A = a 2 kp/2. ) 

stoichiometric constant; ratio of atomic weights of metal to oxygen in the 
oxide (For Cr 2 Og, a = (52.00 x 2)/(16.0 x 3) = 2. 1667.) 

specific linear weight loss rate for the metal in the volatile oxide on the sam- 
pie per unit of time, mg/cm /hr 

stoichiometric constant; ratio of molecular weight of oxide to oxygen in the 
oxide (For CrgOg, b = (52. 00 x 2 + 16. 0 x 3)/(16 x 3) = 3. 1667. It can be 
shown b = a + 1. ) 

conventional parabolic scaling constant in terms of weight of oxygen picked up, 

k„ = — , (mg/cm 2 ) 2 /hr 
P „2 


k y linear weight loss rate due to vaporization of the oxide formed on the sample, 

k = -- x B, mg/cm 2 /hr 
a 

t time, hr 

t time to reach AW in paralinear oxidation, hr 

tg time for AW to return to zero in paralinear oxidation (i.e. , the zero cross- 

over time) 

W specific weight of metal consumed per unit area due to oxide buildup and scale 

vaporization at any time t, mg/cm^ 

W specific weight of metal oxide formed and retained on the metal sample to any 

time t (termed W in ref. 3), mg/cm 

r 

2 

W r lim limiting (maximum) value of W r in paralinear oxidation, mg/cm 
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AW 

AW 

AW 


1»*1 


AW 2 ,t 2 


m 


P r 


o 

specific sample weight chaise at any time t, mg/cm 

maximum AW reached in paralinear oxidation (termed W in ref. 3) 

data sets of observed AW's and their corresponding times 

aggregate average density of the metal(s) consumed in the total oxidation 
process, g/cm 2 

Q 

density of the oxide formed and retained in oxidation, g/cm 


PARALINEAR ANALYSIS 

For conventional diffusion- controlled parabolic oxidation, the basic equations are 

AW=k£/ 2 t 1 / 2 (1) 

W m = ^p 72 * 172 (2) 

and 

W r = bk J/ 2 t 1//2 (3) 

Thus W , the critical attack parameter in oxidation corrosion, is implicit in equa- 
tion (1) assuming, as is almost always the case, that the type of alloy and the composi- 
tion of its scale are known. All pertinent corrosion estimates can therefore be derived 
by simple weight change measurements. 

In paralinear oxidation, two simultaneous processes are occurring: diffusion- 
controlled oxide buildup and loss by such linear processes as scale vaporization and/or 
isothermal scale spalling. An equation analogous to equation (1) can be written for spe- 
cific sample weight change in the approximate form based on the two simultaneous pro- 
cesses, 


AW = kp^t 1 / 2 - kjt 


( 4 ) 


(where k^ is a general linear constant), but no analogous expressions can be written 
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directly for W m or W r for this paralinear model. 

Wajszel, however, in a little noted paper (ref. 3) has started from the basic mass- 
balance equation for any oxidation process. 


AW = W - V/ (5) 

r m ' ' 

and solved the basic differential equations for W r and W m separately. These two 
values can then be combined in equation (5) to give the AW for any set of paired solu- 
tions for W r and W m at any time. Therefore, for the basic paralinear model, this 
method gives the exact solution as a function of time for AW, W m , and W r . 

The basic equations for W m , W p , and AW derived by Wajszel can be generalized 
by dimensional analyses (ref. 9), as described in appendix A. This approach makes 
computer analysis of the paralinear process reasonably straightforward. 

The fundamental differences between the two basic oxidation reactions - parabolic 
and paralinear - are shown in figure 1 in terms of AW, W m , and W r for hypothetical 

alloys with identical k 's but with different k 's. For example, a =2.0 and k = 

9 0 P * p 

0. 1 (mg/cm ) /hr for both; but in one case, k = 0. 05, while in the other case, k = 

p v V 

0. 00001 mg/cm /hr. The latter k y is approximately 0, or the parabolic case in the 
limit. The basic differences are shown on the plot. Not only are the curve shapes 
drastically different, but more importantly in terms of the specific weight of metal con- 
sumed W , in the paralinear case the value at 1000 hours is almost twice that of the 
parabolic case. This difference is even more extreme at longer times. 


FOUR TYPES OF INPUT DATA FOR PARALINEAR OXIDATION SOLUTIONS 

Based on the solution to the equations outlined in appendix A, four cases can be pre- 
sented that cover the major types of analyses for paralinear oxidation data. It is as- 
sumed that the major type of oxide formed has been identified and that therefore the 
stoichiometric constants a and b are known. For the four types of input data de- 
scribed here, all paralinear constants of interest can be calculated - including kp, k y , 
and crossover times. Values for AW, W , and W r can be computed for any desired 
time. The weight change values can be converted to estimated thickness changes when 
the densities of the metal and scale are known. 

The four types of input data are 

Type 1 - The input values are the maximum positive specific weight change AW and the 
time to reach that value t . 

Type 2 - The input values are the parabolic scaling constant k p and linear vaporization 
constant k . 
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Type 3 - The input values are either kp or k y and any one set of weight change and 
time point (e.g. , AWp tj). 

Type 4 - This is the most general case. Two sets of weight change and time points 

(AWptj and AWgjtg) are the necessary input. This case is commonly used 
when a specific we ight- c hange - against- tim e plot is available that extends well 
into the negative weight change region. Here the input can include the last data 
point along with an estimate of the exact time when the AW equaled zero (i. e. , 
the zero crossover point). 


FORTRAN PROGRAM (COREST) 


The program, shown schematically in figure 2, consists of the following 14 sub- 
routines: 

(1) MAIN is the executive routine for processing multiple cases and has primary 
control of logical flow throughout the complete program. 

(2) BLOCK DATA reads initial conditions, parameters, and other data required for 
internal use by the program (looping indices, etc. ). 

(3) BEGIN has charge of initial computations for use by the rest of the program, 
reads program descriptors from cards, and does initial printouts. 

(4) STOICH calculates the a and /3 coefficients for the program. 

(5) ZERO calculates the dimensionless crossover time and the time ratio value of 
the crossover point to the point of maximum net weight. 

The following four QUERY routines determine the numerical values of the paralinear 
parameters from among several possible types of input data: 

(6) QUERY1 - ITYPE = 1; input arguments, AW and t 

(7) QUERY2 - ITYPE = 2; input arguments, k p and k y 

(8) QUERY3 - ITYPE = 3; input arguments, kp or k y and any weight- change- 

against-time pair (AWp tj) 


(9) QUERY4 - ITYPE = 4; input arguments, two weight-time pairs (AWptjj 
AW 2 ,t 2 ) 


(10) OUTPUT calls subroutines for computing the net specific weight curve and con- 
trols printout of results. 

(11) SPEVAL calculates certain weights at special times of interest: W m , W m Q , 

w r , w r lim 

(12) PROVAL calculates the times at which the sample reaches required fractions 
of the limiting specific scale weight. 

(13) CORODE calculates the times at which the sample suffers certain levels of 
metal loss by cdrrosion. 
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(14) CURVE, a principal subprogram, computes the corrosion history of the 
sample. 

The FORTRAN program is listed in appendix B, FORTRAN symbols are defined in 
appendix C, and an example for an experimental case based on type 4 input is presented 
in appendix D. For this example, the total computer time was 23. 6 seconds on a 
UNIVAC 1110. 


USING THE PROGRAM 

Data input to the program is primarily through the BLOCK DATA subprogram 

(fig. 3). After the card "Input follows, " come the actual block data input cards for a 

given series of runs. The contents of each of these cards are as follows: 

Card 1 - DATA ASMALL, BSMALL: The stoichiometric constants a and b for a given 
oxide are listed here. The a value is available in table I for most oxides of 
interest; b = a + 1. 

Card 2 - DATA ITYPE, LL, LLX: The ITYPE values of 1, 2, 3, or 4 specify one of the 
four types of paralinear solutions desired. LL denotes the number of time 
values (up to 100) for which W , AW, W r , etc. , are to be calculated. LLX 
specifies the number of time values (to 10) to be calculated to reach certain 
preselected values of W . 

Card 3 - DATA NUMSET: The number of rims to be computed (to 100) is designated. 

Card 4 - DATA QP: A series (to 10) of decimal fractions of W r are given for 
which the corresponding times are to be computed. 

Card 5 - DATA XXX: A series of values in mg/cm are given for which the correspond- 
ing times in hours to reach these values are to be computed. The number of 
these values listed is equal to the LLX value on card 2. 

Card 6 - DATA HOURS: The individual time values in hours are listed for which the 

corresponding AW, W , W p , and thickness changes are to be calculated. The 
number of values listed is equal to the LL value on card 2. 

Cards 7 and 8 will depend on the DATA ITYPE of card 1. 

Card 7 (type 1) - DATA WBAR: A sequential series of AW values in mg/cm are 

given, with the number of values equal to the NUMSET value on card 3. 

Card 8 (type 1) - DATA TBAR: A sequential series, to match card 7, of "t values in 

hours are given, with the number of values equal to the NUMSET value 
on card 3. 

2 2 

Card 7 (type 2) - DATA XKP: A sequential series of k p values in (mg/cm ) /hr is 

given, with the number of values equal to the NUMSET value on card 3. 
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Card 8 (type 2) - DATA XKV: A sequential series, to match card 7, of as many as 

2 

100 k y values in mg/cm /hr are given, with the number of values 
equal to the NUMSET value on card 3. 

Card 7 (type 3) - DATA NUMKP or DATA NUMKV: Contents are similar to those of 
card 7 (type 2) and card 8 (type 2). 

Card 8-1 (type 3) - DATA WW: A sequential series of single AW values in mg/cm 

at any time are given, to be used with either k p or in the same 
sequence. 

Card 8-2 (type 3) - DATA TT: A sequential series of single-time values t in hours 
are given which match the corresponding AW values on card 8-1. 
These sets of AW, t values are used with the corresponding k p or 
k y . The number of corresponding data sets is equal to the NUMSET 
value on card 3 multiplied by the NUMKP or NUMKV value on card 7. 
Card 7 (type 4) - DATA WW: A series of paired AW values (as many as 100) in 

mg/cm^ at any two times are given, with the number of paired values 
equal to the NUMSET value on card 3. 

Card 8 (type 4) - DATA TT: A series of paired time values t (as many as 100) in hours 
are given for the corresponding AW values on card 7. The number of 
paired values is equal to the NUMSET value on card 3. At program 
execution, five additional input cards are read. 

Card 1A - TRUE or FALSE is punched starting in column 1. FALSE means that the 
term "vaporization" will be used in the headings. TRUE means that the 
term "spalling" will be used in the headings. 

Card 2A - (Format 6A6, F6. 3): A description of the alloy type(s) for the series of 

run(s) printed at the heading of each run is followed by the aggregated aver- 

O 

age density of the alloy being consumed, in g/cm . This density is used to 
calculate the effective thickness change of the alloy based on attack from one 
side. 

Card 3A - (Format 6A6, F6. 3): A description of the scale type(s) for the series of 

run(s) printed at the heading of each run is followed by the estimated density 

q 

of the oxide remaining on the alloy, in g/cm . This density is used to cal- 
culate the scale thickness of the oxide. 

Card 4A - (Format 8A6, 2A6): A description of the data source(s) and the series of 
runs is printed at the heading of each run. 

Card 5A* - (Format 2A6, 2A6): A description of the test temper ature(s) and test pres- 
sure^) is printed at the heading of each run. 

*Card 6A - Plotting routine called (T) or not called (F) with subroutine supplied by 
user at particular installation. 
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An example of this input is shown in figure 3 for a type 4 input of two rims. This 
means there are two sets of paired AWp t^ and AWg, tg input values. There are 17 
time values ranging from 0. 31 to 236. 11 hours for which the AW, W , W r , and thick- 
ness change calculations are to be made. The times to reach three different fractional 
specific weights (0. 5, 0. 90, and 0. 95) of the limiting- scale specific weight are re- 
quested, as well as the times to reach three values of W . The output listing is pre- 
sented in appendix D. 


APPLICATION OF COREST TO THE OXIDATION BEHAVIOR OF 
Cr 2 0 3 /CHR0MITE- FORMING ALLOYS 

The basic COREST program has been used previously by NASA authors with type 3 
input (with a set of AWj,tj values and an estimated k y ) to analyze CrgOg/ chromite- 
forming alloys tested isothermally at 1150° C for 100 hours (ref. 10). It is inferred in 
this reference that under these test conditions the alloys were oxidizing according to the 
paralinear model but that AW occurred well beyond the 100- hour limit of the tests con- 
ducted. Even though the degree of fit appeared excellent, it was desired to examine test 
data where AW had actually been reached and preferably where tQ, the zero crossover 
point, had actually been observed. The literature for classic Cr gOg/chromite- forming 
alloys (i.e. , class n formers in ref. 10) showed two such examples (refs. 4 and 6). 

For this reason, four supplemental tests were performed on two additional alloys at 
this laboratory for times long enough so that at least the zero crossover point was ap- 
proached. This testing was conducted basically to determine how well the early (<t ), 
middle (~T), and later (>t) time parts of the curve establish the rest of the AW- 
against-time plot and how their use changes the long-time predictions (>10 000 hr). 

The paralinear model, assuming it still holds, is thought to be more realistic than the 
parabolic model. 

The data from the two literature and four NASA Lewis runs used for input into 
COREST are shown in table II. Data on Cr-0. 2Y were taken from reference 4, and data 
on Ni-20Cr-3Fe-0. 4Si (termed ’’Alloy A") were taken from reference 6. Both alloys 
were tested in static air at 1200° C. The tests performed at this laboratory were on 
Ni-40Cr and Tophet 30 (Ni-29. 6Cr-l. 4Si) at 1200° C in static air. In addition, duplicate 
tests were run on Tophet 30 at 1150° C in static air. 

The six experimental AW- against- time plots are shown in figures 4 to 9 with the 
observed data points plotted as circular symbols. Also shown on each plot (square sym- 
bols) are the COREST-generated AW- against- time plots derived from the type of input 
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data that gives the lowest error of estimate. 

Table HI shows the degree of fit based on the minimum error of estimate calculated 
between the observed values and each comparable COREST-calculated data point. Also 
listed are the various paralinear parameters of interest: kp, k y , AW, and tQ and the 
predicted t for W m to reach 71. 9 mg/cm 2 (i. e. , 100 pm of metal consumption). 
Finally, the critical W m values are plotted in figure 10 for times longer than 
1000 hours for each of the six alloys. The W m values plotted for each alloy are those 
representing the best degree of fit. 


RESULTS AND DISCUSSION 


The specific weight change curves for the six experimental runs are plotted in fig- 
ures 4 to 9 as a function of time, along with the COREST-calculated AW values. These 
runs are summarized in table in, which lists the three time stages of input data from 
table n, as well as the various derived paralinear constants and other values of interest 

such as selected estimates of certain t and W„ values. The degree of fit between the 

r 2 

derived and actual curves is given as an average error of estimate in mg/cm . These 
cases were derived either from a type 1 calculation, where the curve maximum was 
estimated, or from the more general type 4 calculation used both for the pre- and post- 

O 

AW cases. The degree of fit is good - usually less than ±0. 5 mg/cm even for the 
longer runs. 

The degree-of-fit criterion is given in terms of an average error of estimate of the 
calculated AW curve compared with the experimental AW curve. It was used to deter- 
mine how long, in general, a paralinear test must be run to give meaningful long-time 
estimates of W . Obviously, from a practical point of view, a shorter (before AW) 
100- to 200-hour test is to be preferred. However, as the degree-of-fit criterion indi- 
cates, the tests should at least be run to AW - and preferably into the negative slope 
region. Using input values before AW does, however, at least result in more pessi- 
mistic long-time W m estimates than the AW or post- AW cases. On the other hand, 
if parabolic oxidation is assumed, which has usually been the case, the long-time W m 
values are drastically underestimated, as shown in figure 10. 

An example of how the assumption of parabolic kinetics can influence the interpreta- 
tion of paralinear oxidation data is shown in figure 11. If the Alloy A oxidation test dis- 
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E 


< AW i,exp- Aw i,calc> ‘ 
Number of AW values 


9 * 11/2 


= Average error of estimate 
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cussed previously had been only 300 hours long rather than the actual 1000 hours, the 
resulting data plotted on a (AW) -against-time plot would appear to represent two-stage 
parabolic oxidation. The derived parabolic constants, if then used to estimate metal 

o 

loss, would give a W m at 1000 hours of about 16. 5 mg/cm . If only the higher first- 
stage constant were used, metal loss would be about 20.7 mg/cm 2 . However, based on 
the full 1000- hour curve with the vaporization of Cr 2 Og leading to paralinear behavior, 
the W m estimate is 32.61 mg/cm 2 , almost twice the original estimate. Extrapolated 
to 10 000 hours, this ratio of underestimated loss would be almost 5:1. On the other 
hand, if short-time data were used in a type 4 paralinear analysis, a greater W m value 
would actually be estimated at 1000 hours than when the negative slope values were used. 
This is the most conservative estimate and probably the safest from a design standpoint 
if the W m values are converted to an effective thickness change. It has the obvious 
advantage that a test would not have to be run out to a weight loss situation. 

The long-time comparisons for Alloy A using all three of the basic attack vari- 
ables - AW, W r , and W m - are shown in figure 12 for both the parabolic and paralinear 
(calculated) cases, assuming the same kp. This log- log plot graphically shows that as 
the scale specific weight approaches about 90 percent of its limiting scale specific 
weight of 64. 5 mg/cm , the W m value rises sharply and the AW drops sharply (to- 
ward negative values). 

The calculated kp values (based on input values with the best fit) are replotted as 
a function of the reciprocal of the absolute temperature in figure 13. The kp values 
fall in the range for conventional CrgOg-forming alloys, which is discussed in refer- 
ence 10. These were also calculated by COREST but based mostly on shorter tinle 
(100 hr) exposure data. The kp results for Tophet 30 fall slightly below the lower limit 
for the Cr gOg/chromite range. It is believed that the relatively high Si content in 
Tophet 30 causes the erratic early paralinear oxidation behavior as well as the lower 
consumption rate after AW. Silicon dioxide (SiC> 2 ) in the scale apparently behaves like 
thorium dioxide (ThC> 2 ) blockage in TD-NiCr alloys, lowering the effective kp (ref. 10). 

The derived paralinear k y values (again based on input values with the best fit) are 
plotted in figure 14 against 1/T K , along with other pertinent k y values derived for the 
vaporization of bulk Cr 2 Og. This figure shows k y values derived gravimetrically by 
Hagel (ref. 11) for sintered Cr 2 Og tested in still air, along with the limits for 1 stand- 
ard deviation. The k y value of a sintered Cr 2 Og sample run at this laboratory as a 
check falls almost exactly on Hagel’s value. This implies that test conditions were com- 
parable and that condensation of Cr 2 Og on the balance hanger wire did not significantly 
bias the derived results. Also plotted in figure 14 are literature values for vaporization 
of Cr 2 Og prepared by sintering or hot pressing both in static and moving air (ref. 12) 

o o 

The slope of the (AW) -against-time plot is the usual graphical method used to de- 
rive parabolic scaling constants. 
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and for vaporization of stripped Cr 2 Og scales from an oxidized Ni-50Cr alloy (ref. 8). 
Ecer and Meiser (ref. 8) also calculated a minimum vaporization rate in static oxygen 
at 1 atmosphere based on gaseous diffusivity. Of interest from a practical point of view 
is the maximum possible loss rate of Cr^Og. Such loss can be determined from the 
vapor pressure of CrOg at the partial pressure of oxygen in 1-atmosphere air by using 
the Langmuir equation and vapor pressure data for CrOg cited by Kohl and Stearns 
(ref. 13). Obviously, this is a limiting case with a Cr loss rate some 4 orders of mag- 
nitude higher than that observed in static air. The ky loss rate in a high-velocity 
burner rig test on a CrgOg-forming alloy should be intermediate between these limiting 
values but as yet is undetermined. As the upper limit is approached, the Cr consump- 
tion rates would be so high, particularly in the higher temperature regions, that the Cr 
would probably be depleted, which would favor NiO and increased chromite spinel forma- 
tion. In cyclic oxidation these scales would spall catastrophically (ref. 14), which would 
result in accelerated attack. 

It can be seen that these calculated ky values can vary widely at a given tempera- 
ture between the upper and lower estimated limits - over 5 orders of magnitude. Most 
laboratory tests, whether in still or moving air, appear to be within a factor of 10 in the 
lower third of the range. It was thought that, in the static-air tests, the paralinear k y 
values were lower than Hagers bulk values because NiO and/or spinel overlay the Cr 2 Og 
(ref. 10). However, the scatter may be due to experimental variation. It does appear 
that the derived ky values for the three Tophet 30 runs are somewhat lower than those 
for the other alloys cited. If these are real differences, the implication in terms of 
overall attack (i. e. , W m ) is that Tophet 30 is even more oxidation resistant in static air 
than the general class of Cr 2 Og/chromite-forming alloys. 

It is critical, particularly for a practical test situation, to determine how long a 
test exposure should be run before deciding that the oxidation behavior is paralinear 

rather than parabolic. For Cr 0 0 Q /chromite-forming alloys, figure 13 shows that k„ 

^ o _ o q no P 

values can range from just over 10 to just over 10 (mg/cm ) /hr, in the temperature 
range 1000° to 1200° C. From figure 14 the lower range of ky values in static air is 
from about 10" ^ mg/cm^/hr at 1000° C to as high as 1. 0 mg/cm^/hr at 1200° C. These 
ranges of kp and k y values were combined in a series of COREST calculations to pre- 
dict "t, the time to reach AW. This is the minimum test time to determine paralinear 

behavior. Because the inflection tends to be gradual and AW returns to the origin 
o 

(0. 0 mg/cm ) at just over four times the predicted t, running the test to twice the pre- 
dicted T seems reasonable to accurately determine the behavior. 

This approach is summarized in figure 15, which shows the critical nature of k . 

O o 2 ” 

values between 10 and 10 mg/cm /hr. The series of curves imply that for virtually 

_3 

all k values of 10 or less the behavior is effectively parabolic for most practical 
v -2 2 

exposure situations. Above k y values of about 5x10" mg/cm the kinetics are 
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obviously paralinear for test times of just a few hundred hours. For k values between 
10” and 5x10 mg/cm /hr, figure 15 can be used in conjunction with the end applica- 
tion to determine a meaningful test time. Obviously, for tests in a moving- air environ- 
ment the k y values shift to the higher values in figure 14, where paralinear behavior is 
apparent at much shorter test times. 

Figure 16 shows how this paralinear approach might be used to estimate long-time 
oxidation damage over a wide range of temperatures. From figure 13, upper and lover 
limit values of kp for predominantly Cr 2 Og-forming alloys were estimated at 50- 
degree intervals from 1000° to 1200° C. From figure 14, for each of these tempera- 
tures three values of k were estimated, representing the likeliest static-air value, 

10 times this value (representing moving air), and the maximum possible k y derived 
from the Langmuir equation for free vaporization (ref. 13). These values of kp and k y 
were combined in the six possible combinations as the input values in the COREST pro- 
gram to generate AW, W m , and W r as a function of time. A plotting subroutine was 
then used to plot W m as a function of time on a log- log plot at times longer than 
10 000 hours. The pertinent data (except for the free- vaporization case) are summar- 
ized in table IV. These curves are based on the assumption that paralinear behavior 
continues beyond 10 000 hours, which implies continued Cr 2 Og/ chromite control. How- 
ever, if the base metal oxide becomes rate controlling because of previous Cr depletion, 
earlier catastrophic failure is apt to occur from spalling - particularly if thermal 
cycling is involved. The paralinear estimate is still more realistic than that derived 
from the parabolic model. However, comparable parabolic plots on a log-log scale 
would be straight lines with a slope of 0. 5. 

At all temperatures, when the Langmuir vaporization rate is the controlling factor, 
the overall loss rate is effectively a linear process. In this case the ratio of k y to kp 
is so great that the scale has a limiting thickness of well under a micrometer. This 
ratio is close to 0. 01 with the low- range k^'s and equal to 0.25 with the high- range 
kp's. When the ky value is increased by a factor of 10, the limiting scale thickness 
drops by a factor of 10, and the time to consume a given amount of Cr decreases by a 
factor of 100. Thus, the metal consumption rate seems to increase with the square of 
the ratio of the vaporization rates. Parabolic scaling constant ratios, on the other hand, 
affect metal consumption rates in strictly a linear manner. The times to reach a given 
W_, show the striking differences at all temperatures between the paralinear and para- 
bolic cases. Ratios are as high as 39:1 for W m =719 mg/cm (~ 1000 jum) for even the 
static-air case. 


CONCLUDING REMARKS 

The oxidation behavior of four predominantly Cr 2 Og-forming alloys oxidized for 
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long times above 1100° C was analyzed by comparing the observed net- weight- change - 
against-time values with their corresponding calculated values, assuming paralinear 
behavior. The calculated values were generated by the COREST computer program 
based on short-time input before the obvious AW, at the estimated AW, and after the 
AW on the downslope of the curve, where paralinear behavior is clearly established. 

The agreement between the calculated and observed AW values on the average was 

O 

±0. 5 mg/cm . Based on the degree of fit as might be expected, long-time input data are 
preferable. However, even the use of short-time input data, assuming paralinear be- 
havior, is more meaningful than assuming parabolic scaling kinetics. The latter method 
greatly underestimates the critical W m values particularly at longer times. Based on 
these results, it is suggested that the paralinear assumption be used to analyze any 
short-time data since it includes the parabolic case as a limit with k„ — 0. The k„ is 

p V V 

effectively zero if it is <0. 001 mg/cm /hr, provided k„ is greater than about 
2 2 P 

0. 01 (mg/cm ) /hr. In this case the behavior is essentially parabolic and the test can 

be terminated. 

Based on the range of derived kp and k y values for CrgOg/ chromite-forming 
alloys, a series of curves have been derived for the critical corrosion parameter W m 
at 1000°, 1050°, 1100°, 1150°, and 1200° C for times to 10 000 hours and beyond. The 
six curves generated at each temperature represent the highest and lowest kp values, 
along with three k y values - the likeliest static-air value, 10 times this value (which 
represents moving air), and a maximum k y derived from the Langmuir equation for 
free vaporization. The curves reflect the strong effect of vaporization, particularly in 
moving air. 

The COREST paralinear mass-balance approach is currently being used in analyzing 
cyclic oxidation where scale spalling, as well as oxide vaporization, can occur. The 
COREST approach appears quite promising for those alloys where the spalling rate ap- 
pears fairly uniform over a number of test cycles. 


SUMMARY OF RESULTS 

A FORTRAN computer program termed COREST has been developed for analyzing 
paralinear oxidation based on the mass-balance expression AW (net sample weight 
change) = W y - W . The W y (specific weight of the retained oxide) and W m (specific 
weight of the metal consumed) values can be generated independently as a function of 
time from a minimal amount of input data with a known scale. The W m value is 
critical in corrosion, and the W m values in paralinear oxidation are always increas- 
ingly greater with time than the corresponding values for parabolic oxidation. The 
COREST program also converts W y and W m to thickness values with suitable input of 
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oxide and metal density values. With a stoichiometric constant calculated from the com- 
position of the oxide, the following four types of input data can be completely analyzed 
by COREST: 

(1) The maximum estimate value of AW and its corresponding time 

(2) Any paired values of the parabolic scaling constant kp and linear vaporization 

rate constant k y (where k v + 0) 

(3) Any value of kp or ky, along with any one set of AW with its corresponding 

time 

(4) Any two sets of AW values with their corresponding times 
From these four types of input the program calculates and lists 

(1) kp and ky 

(2) The maximum positive AW and its corresponding time and W m and W r 

(3) The time for AW to return to zero and its corresponding W m and W r 

(4) The limiting value of W f and its corresponding thickness 

(5) The time(s) to reach any requested fraction(s) of the limiting value of W r 

(6) The time(s) to reach any requested value(s) of W m 

(7) The AW, W m , W r , effective metal thickness loss, and oxide thickness for any 

requested time(s) 

The derived paralinear scaling constants for straightforward chromic oxide 
(CrgOg )/ chromite-forming alloys generally fall within the range of the previously pub- 
lished values of kp. The derived k v values increase as expected with increasing tem- 
perature but are less than the k y values for bulk CrgOg in static air. This difference 
is greater in NiCr alloys than in pure Cr presumably because of a thin NiO and/or 
chromite spinel layer on the outer surface which probably spalls on cooling. 

The presence of 1. 5 percent silicon in Tophet 30 appears to lead to silicon dioxide 
blockage in the CrgOg scale. This apparent blockage causes the initial kp to drop until 
the time to reach AW, the maximum AW reached in paralinear oxidation, is exceeded 
and thus leads to lower than expected metal consumption rates. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 30, 1975, 

505-01. 
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APPENDIX A 


ANALYTICAL EQUATIONS AND COMPUTER PROGRAM 


Derivation and Solution of the Paralinear Equations 


The analysis anticipates a type of metal oxidation experiment from which the basic 
data consist of a history of sample weights taken at a progression of discrete times and 
a knowledge of the stoichiometric conversion factor from metal to oxide scale. For 
those systems that suffer scale loss by vaporization, the mass balance for the sample 
states that 


'Total sample 
.weight change. 


Weight of 
.scale produced. 


Weight of 
.metal oxidized. 


Weight of 
.scale vaporized. 


(Al) 


Now, for the metal- oxidation reaction, 


mM + ( f )°2 


M^ C> 
m n 


(A2) 


the stoichiometric conversion factor a is defined as 


a = - 


Atomic weight of" 
m 


metal, M 


Molecular weight of 

- m 

Atomic weight of 

scale, M m O„ 

L ’ m n J 

metal, M 


(A3) 


As an example, a = 2. 1667 for chromic oxide (CrgOg). The ratio of the weight of scale 
formed to the weight of metal oxidized is 


"Molecular weight of"! 

scale, M 0„ 

* m n J 

"Atomic weight of] 
m 

metal, M J 


(A4) 


Replacing the scale weight term in the mass-balance equation (Al) by the oxidized metal 
weight times the ratio (A4) yields 
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' Total sample ' 
weight change. 


g 1 

B 

1 

^ i 

r Weight of 1 _ [" Weight of 

L mM J 

Lmetal oxidizedJ Lscale vaporized. 


1 1" Weight of 
a Lmetal oxidized. 


_s 


Weight of 
cale vaporize 


d_ 


The time derivative of equation (A5) 


dW _ 1 dW m 
dt a dt 


dW^ 

dt 


(A5) 


(A6) 


is the mass-balance rate equation that is the starting point for subsequent analyses of 
the metal oxidation history. 

By using the parabolic rate term of Haycock (ref. 4) for dW m /dt and a linear 

vaporization rate for dW„/dt, Wajszel (ref. 3) was able to integrate equation (A6). 

s 

This solution for sample weight as a function of time was termed paralinear by Wajszel 
because of its dependency on the parabolic oxidation model and the linear scale vapori- 
zation model. 

Paralinea r oxidation model . - This is an Ohm's law analogy. The flux of oxygen 
ions through a unit cross-sectional area to the scale-metal interface is inversely pro- 
portional to the length of the flux path (i.e. , scale thickness). The proportionality con- 
stant is the paralinear conductance of Haycock. 

Since the rate of metal oxidation equals the flux of oxygen ions, Wajszel expresses 
this fact within the Ohm's law model as 


dw m /b\ A 
dt (a/ W r 

4 

where the weight of scale per unit area W r is equivalent to the scale thickness, the 
proportionality constant A contains the paralinear conductance, and the factor b/a is 
the molecular weight ratio of scale weight to metal weight. 

Linear sca le vaporization model . - The simplest model one can devise for loss of 
scale by vaporization postulates a rate that is independent of time for a given set of ex- 
perimental conditions. In other words, the total metal loss through scale vaporization 
is simply linear in time, or 


4 


In 


Wajszel's paper (ref. 3) this is termed W n . 

r 
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= Bt 


(A8) 


^ “'vaporization loss 

where B in the vaporization rate constant and t is elapsed time. Now the scale weight 
loss is simply the weight ratio factor b/a times the left side of equation (A8), or 

W g = ^Bt (A9) 

Integration of pr oduct scale equation. - One form of the mass-balance equation is 

r eight ofl _ b f" Weight of 
scale J a Lmetal oxidize 


d. 


Weight of 
.scale vaporized. 


(A10) 


or 


W_ 



The time derivative of equation (All) is 


dW r 

dt 


h dW 
b m 

a dt 


dW^ 

dt 


(All) 


(A12) 


Substituting equations (A7) and (A9) into the right side of equation (A12) reproduces 
Wajszel's equation (5), 


dW r 

dt 


b 

a 



(A13) 


It is seen from equation (A13) that when the derivative vanishes, we have the maximum 
(steady state) scale formation, that is, 

(W ) = ^ (A14) 

'max a B 


which will be shown to occur from the integrated form of equation (A13) as time ap- 
proaches infinity as a limit. 
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Equation (A13) takes the form 


— g dz =dt 

B 2 (A - z) 

where z = BaW r /b. Making a second change in the variable yields 

u = A - z 
z = A - u 
dz = -du 

Then 


(A15) 



This finally is rearranged to Wajszel's equation (8) 



(A 16) 


We find the exponential form of equation (A16) to be more useful in the subsequent 
analysis: 


R = 1 


-R-T 


(A17) 

19 



where 


aBW_ 

R = ' 

bA 


and 



(A 18) 


Integration of the metal- oxidation equation . - When the scale product equation (All) 
is written with equation (A9) substituted for the vaporization term, the resulting equation 



- Bt) 


(A19) 


can at once be substituted into the right side of equation (A7) to obtain the important 
metal- oxidation equation 


dW. 


m 


dt W. 


m 


Bt 


(A20) 


The last expression is Wajszel's equation (3a), which, along with equation (A13) (or 
Wajszel’s (5)), constitutes the analytical basis for the computation program of this re- 
port. 

The integration of equation (A20) will now be demonstrated. First, make the trans- 
formation of the variables 


M = 


W m B 


T = 


tB 


J 


whereupon the new form of the differential equation becomes 


dM 1 
dT M - T 


(A21) 


(A22) 
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Further, let M - T = 1/v. The limits on these variables are 


for T=oJ< M - T >=° 

[v — °° 

for T - » (< M - T > - 1 
It- 1 


The integral of equation (A22) takes the form 


J-oo V 2 (V 2 - 1) J 0 


dT 


Let the integrand on the left be expanded into the partial fraction 


-1 = 1 + 1 1 


2/ - v v 2v-l 

V (v - 1) V 


This allows equation (A24) to be immediately integrated 


fin v - v _ * - In (v - 1) = T 

L J-oo 


or 


In 


v 1 


V - 1 V 


= T 


And since the bracket on the left vanishes at the lower limit, we then arrive at 


In — - -i = T 


v - 1 v 


or 


In 


1 - M + T 


- (M - T) = T 


(A23) 


(A24) 


(A25) 


(A26) 


(A27) 


(A28) 


(A29) 


21 



and finally 


In (1 - M + T) = -M 


(A30) 


Equation (A30) is one form of the definite integral of equation (A20) in terms of the non- 
dimensional variables M and T. 

The exponential form of equation (A30) is more useful in the computer program. So 


M = 1 + T - exp (-M) 


(A31) 


which is exactly Wajszel's integral solution (eq. (7)) in terms of M and T. Replacing 
the M and T variables in equation (A31) by their definitions (A21) yields 


B 


-1 + exp (- 


■?) 


W m 
m . 
+ = t 

B 


(A3 2) 


which is Wajszel’s equation (7). 

Estimation of sample weight change . - It is now possible to derive an expression for 
calculating the predicted specific sample weight change due to oxidation of the metal. In 
the mass-balance equation (Al) the term for the total scale formed minus the term for 
scale vaporization is just the product scale weight. Thus, the mass equation may be 
written as 


'Specific sample - ] [Product scaled T Weight of 
_ weight change J L weight Lmetal oxidized. 


(A3 3) 


or 


AW = W r - W m (A3 4) 

The terms on the right are obtained from equations (A16) and (A32), respectively. 

When numerical values are given for the system parameters a, A, and B, the 
transcendental equations (A16) and (A32) can be solved for W m and W p at any time t. 
It is necessary to solve equations (A16) and (A32) iteratively. The calculations to 
achieve this are described in the next section. 
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Numerical Procedures for the Paralinear Equations 


Dimensionless for ms of the equations . - The integral solutions of the W m and W r 
variables have already been obtained in equations (A16) and (A32). The transformed, 
nondimensionalized variables introduced in definitions (A18) and (A21) then restated the 
solution equations for W m and W r in the nondimensional forms (A17) and (A31). 

The fundamental mass-balance equation (A34) will be transformed in like manner. 
Multiply by the ratio of coefficients B/A 


B AW 
A 


BW r 

A 


BW 
m 

A 


(A3 5a) 


B AW b /a B \ BW m 


which is reduced to 


N =-R - M 
a 


(A35b) 


(A3 6) 


In this equation the transformed, nondimensionalized specific sample weight change N 
is the dependent variable as a function of the scaled independent variable time T. 

The phrasing of physical laws in terms of nondimensional variables, as far as pos- 
sible, is a rather common practice within disciplines which use modeling techniques to 
cope with very complicated problems. Among these other fields are aerodynamics and 
hydraulics, hydrology, chemical reactor design, river and harbor studies, and analysis 
of smoke dispersion from stacks. The technique is known by a variety of labels: simil- 
itude, modeling and scaling, similarity analysis. In its crudest form, it is represented 
by the dimensional analysis used by Lord Rayleigh and Osborne Reynolds (ref. 9). 

The present study has already demonstrated that a practical advantage of scaling 
(nondimensionalizing) arises from the simplification of the form of the model equations. 
It will be apparent later in the discussion of the iterative solution of equations (A17) and 
(A31) that the computer subprogram requires no transfer of system coefficients or 
parameters other than the independent variable T. Also, further insight into the impli- 
cations of the basic system models results from working with the simpler forms and 
solutions of the scaling analysis. This will now lead to some additional and useful prop- 
erties. 

Properties an d trends of t he model . - The mathematical limits of the relations 
(A17), (A31), and (A36) are 
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fR = 0 
for T = 0<M = 0 

In = 0 


(A37) 


for T >> 



Equation (A36) will be examined for an optimum point. But first note that equations (A17) 
and (A31) have equivalent forms if 


M = R + T (A3 8) 

which is actually a restatement of the mass-balance equation (A10). To search for a 
maximum point on the (N, T) curve, the derivative of equation (A36) is set to zero, 


Then 


and 


dN _ b dR dM _ Q 

(A3 9a) 

dT a dT dT 


bl-R. 1 =0 

a R M - T 

(A39b) 


bl-R 1 _ . 

(A39c) 


a R 


R 


The terms for dR/dT and dM/dT on the right in equation (A3 9b) are from equations 
(A17) and (A31) and use has been made of equation (A38) to obtain (A39c). The special 
notation for R, M, and T refers to those values at optimum N. 

Equation (A39c) results in 


b 

The logarithmic form of equation (A17) shows time at maximum N to be 

T = -R - In (1 - R) 



(A40) 


(A41) 
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From equation (A3 8), 



and for the maximum N, from equation (A3 6), 



(A 42) 


(A43) 


It is of interest to note that the values of N, M, R, and T depend only on the value of 
the stoichiometric coefficient a. These special values for N and T have been given 
the notation 


a 



'--OH J 


(A44) 


From the definitions of T and N in equations (A18) and (A36), respectively, it 
can be shown that 


and 



(A45a) 


(A45b) 
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from which it should be noted that the scaling of N and T is based on AW and t, 
respectively, the maximum specific net sample weight change and the time at which it 
occurs. 

The rate coefficients A and B reflect the values of AW and T. From the defi- 
nitions of N and T (eqs. (A36) and (A18)) and the optimum values N = a, T = /3, there 
results 


A = 


B = 



(A46) 


This last relation provides a most convenient method for calculating the rate coefficients 
from limited experimental data. 

Another obvious examination point occurs at the time Tq when N = 0, which is 
designated the ’’crossover point. ” From equation (A36) we have 


— - = M 0 (A47) 

a 

Using equation (A3 8), this becomes 

Rq = aTQ (A48a) 

or 

M Q = bT Q (A 48b) 

Either of the forms of equation (A48) substituted into equation (A17) for Rq, or into 
equation (A31) for Mq leads to the equation 


aT 0 = 1 - 


-bT 

e 


0 


(A 49) 


which must be solved by an approximate technique. This also shows that Tq is only a 
function of the stoichiometric coefficient. Table V illustrates the general values of 
these dimensionless variables over a range of values of the coefficient a. 
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One of the consequences of the paralinear model is obtained from the last column in 
table V, which shows that the ratio of crossover time to time at optimum sample weight 
change Tq/T is about 4. 

Table VI summarizes the results of this analytical section on the properties and 
trends of the model. Four time values are exhibited: initial time t = 0, time at maxi- 
mum specimen weight gain t, time at zero specimen weight gain t^, and finally very 
long times where the variables N, M, R, and dN/dT approach limiting values. 

Table VII summarizes the definitions of the much-used paralinear parameters and 
their development from basic physical variables (i.e. , their correspondence with the 
transport coefficients kp and k y ). 

N umerical solutions of the equations . - We have seen that the nature of the para- 
linear model is such that the history (or curve) of AW as a function of t is uniquely 
predicted from a knowledge of three constants: a, A, and B. 

Conversely, two data points (AWj,t^) and (AW 2 ,t 2 ) can, in principle, be used to 
calculate the values of A and B that are needed to give the unique curve passing 
through those data points. This is, of course, about the only way to obtain estimates of 
A and B from the laboratory experiments. 

This section outlines those numerical techniques for solving the paralinear equations 
that are the algorithms for the FORTRAN program. The complete program procedure 
is in the following sequence : 

(1) Establish values of the system parameters a, y8. A, B, Tq 

(2) Solve the transcendental paralinear equations in dimensionless form by approxi- 

mation methods 

(3) Use these solutions to obtain a group of results: 

(a) The array of AW's for a given array of t values 

(b) The special results W m , Q , W m , W r , W r lim 

(c) The times tg at which W r = 6 • W f ^ 

(d) The times t corresponding to chosen values of W m 

Syste m parameters . - The crossover point Tq is obtained from equation (A49) by 
using successive substitutions of improved approximations to Tq, starting with an ini- 
tial value of Tq = 4/3. The convergence was rapid in satisfying the test 


T, 


j+1 


T. 

J 


< 0.0001 


T i 


(A50) 


The values of A and B are obtained by using one of four data input types: (1) de- 
fined values of k and k (table VI); (2) given values of (AW, t ) (eq. (A46); (3) a given 

P V 

value of kp or ky and one general data pair (AW, t); (4) experimental data pairs 
(AWpt^) and (AW 2 ,t 2 ). The computation using the third type differs depending on 
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I 


whether kp or ky is defined. 
Example 1 - Define k : 

hr 


A = 



2 


From equations (A3 6) and (A3 8), 


R = a(N + T) 

= — (A W+Bt ) (A51) 

A 

The terms in (A51) are known except for B, which makes the equation quadratic in B. 
Both equations (A51) and (A17) are to be satisfied. Write equation (A17) in its alternate 
form 


T = In (1 - R) - R 


(A17a) 


or 


f(B) = T + R + ln(l-R) 


(A17b) 


where f(B) = 0 when the correct value of B is obtained to satisfy both equations. The 
B-root of equation (A17b) is obtained by using the Newton-Raphson technique 


where 



(A52) 


(A53) 
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II 


where 


/dT\ 2B j t 

\ dB j/ A 


— =— (R; + aTJ 
d Bj Bj 1 J 


R = — 1 (AW + Bit) 
J A ] 


B?t 

T = J 
J A 


f(Bj) = Tj + Rj + In (1 - Rj) 

The initial value of B^ is obtained from the definition of T, that is, 



(A54) 


where is arbitrarily taken at Rj = 0. 5, that is, half of its maximum value. Con- 
vergence of the Newton- Raphson method to the limit test is rapid: 


l B i+i - B i i 


< 0. 00001 


(A55) 


Example 2 - Define k y : 



(A46) 


Equations (A51) and (A17b) are again used, but with A as the unknown coefficient. By 
analogy to example 1, 


f(Aj) = T. + Rj + In (1 - R.) 


(A56) 
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where f(A) = 0 at convergence. Again Newton-Raphson uses 


and 


f(A,) 

A, j = A, L 

3+1 3 f'(A.) 

J 


(A57) 


where 


f’( A ) 

3 dA. 


dT 


R j dR 

dA. 1 - R. dA. 

J J J 


dT 

dA j 


T i 


dR 

dA j 


h 

A i 


Rj 

-T, + — L 


f, ( A j) = 


3 1-R. 


A i 


(A58) 


Numerical iteration is handled as in example 1. 

Example 3 - Define two experimental points, (AWptj) and (AWg, tg): We start with the 
mass-balance equation (A36) and replace R by equation (A3 8) 


aN = M - bT (A59a) 

— AW = M - bT (A59b) 

A 
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For both pairs, when the weights are nonzero, we have for the two times (tg > tj), from 
equation (A59b), 


Mj - bTj M 2 - bTg 
AW X ~ AW 2 


/AWA / AW, \ 

M-, - [ i]M 9 + bT 1 |T I =f(T,) (A60) 

V AW 2/ V AW, ) 

where r = Tg/T ^ = tg/tj and f(Tj) = 0 when the parameters A and B are evaluated 
such that and Tg = tT^ give and Mg values which satisfy equation (A60). 

The root of equation (A60) is found by using the Wegstein secant modification of 
the iteration technique for successive estimates on T^: 

T i+ i = «V 


AWjMg 

AWo 


- M, 


' AWjtg 
^Wgtj 


- 1 


(A61) 


where Mg and M^ are calculated for given rT- and T-, respectively, by using the 
Newton- Raphson iterative method to solve equation (A31): 


M. 


i+1 - M i " 


F(M i ) 

F’(M.) 


(A62) 


where 


F(Mj) = -Mj + 1 + Tj - exp (-M^ 


(A63) 


F’(M-) =— = T i - M i 

1 dMj J 1 


(A64) 


The Newton- Raphson method rapidly converges for given values of T-. 

J 
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A safe starting value for Tj in equation (A61) is Tq, from (A50), and the initial 
M 0 = bT Q in equations (A62) and (A64). 

Wegstein’s iteration scheme is used to get accelerated convergence in the case of a 
function, such as equation (A61), that has at least a continuous first derivative in the 
range in which the iteration moves. It refines an initial guess T Q of the root of our 
nonlinear equation T = f(T). 

The calculation sequence starts as follows 

x i = y 0 = f ( x 0 > 

yj = f(x x ) 


with the succeeding steps resulting in the computational scheme 


x i+ l = x. + 


x i - x i_i 


* 1-1 - * 1-1 
x i“ Yi 


- 1 


> a = 1 , 2 ,...) 


yi + i= f(3 4 + i> 


(A65) 


Each step requires one evaluation of f(x). 

The iterative procedure is terminated when the following two conditions are satis 

fied: 


DELI < e and DEL2 < 10 • e 


with 




DELI =< 


*1+1 -*i 


when Ix^j] > 1 


*i+l 

lh+1 ' *il when l x i+ll - 1 
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DEL2 =< 


*4+1 ~ ^i+1 

T i + 1 


when |x i+1 | > 1 


h+i - y i+ ii when i x i+ ii s 1 


and the tolerance factor e is specified in the input. 

The calculation of the parameters A and B follows after convergence of equa- 
tions (A61) and (A64) to a suitable test, such as was used for equation (A54). The com- 
puted values of Mj and Tj are used in equation (A59) to determine 


B M l- bT l 
A a AW 1 


and from equation (A18) 




(A66) 


(A67) 


Dividing equation (A67) by (A66) gives the sublimation parameter 


T i / M i- bT i 

' ‘i V a iW i / 


(A68) 


and finally the parabolic parameter A is obtained from either equation (A66) or (A67), 
for instance, 


A 



(A69) 


For the special case in which one given weight is zero (crossover point), the associ- 
ated time is tQ, which establishes the important T and thus the B^/A = 0/F. The non- 
zero weight-against-time pair (AW, t) is used to establish the time T = (/3/t)t and the 
associated M-value from equation (A64). Finally, from equations (A68) and (A69) the 
individual parameters A and B are computed. 

Table vm summarizes the four possible input types used to establish the necessary 
system parameters. 
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Main program calculation. - With numerical values of the paralinear parameters 
a, A, and B, the main results can be calculated by using the relations already estab- 
lished. 

Net specimen weight curve: Given an array of LL times at which values of net 
specimen weight are wanted, specify t t 2 , . . - , tj^ij hours and calculate 

TpTg,. • . ,T t T from equation (A21) 

• . ,M T T from equation (A3 1) 

Nj,N 2 ,. • . , N T T from equation (A59a) 

AWpAWj,. . . AW^ l from equation (A5 9b) 


W m, 1* W m, 2 ’ ' • •’ W m, LL from equation (A21) 
W r lf W r 2> * • • 5 w r ll from equation (A18) 

6 m, V 6 m, 2> • * * ’ 6 m, LL irom e ^ tion ( A7 °) 

5 r,l’ 6 r,2’* • - } 6 r,LL 

The last two rows of thicknesses are calculated as 


(Metal thickness lost) = — — - 6 (A70) 

p _ 
m 

where metal density is p , and for the thickness of the scale layer at any time 

W„ 

(Scale thickness) = = 6 (A71) 

^r 

where P r is the scale density. These calculations are found in FORTRAN subroutine 
CURVE. 

Special weight values: Table IX completes the trends and limiting values of table VI. 
These calculations are found in the FORTRAN subroutine SPEVAL. 

Product scale thickness-against-time values: It is required to determine the exper- 
imental times for the model sample to reach certain specified fractions of the equilib- 
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rium (limiting) scale thickness. Given the fractions Rp R 2 , R3 < 1, calculate Tp 
T 2 , and Tg from equation (A16) and tp t 2 , and tg from equation (A21). These are 
calculated from the FORTRAN subroutine PROVAL. 

Metal corrosion- time values: Analogous to the preceding paragraph, estimate the 
times required for the sample to suffer corrosion losses of specified values. Thus, 
given LLX values (W m p W m 2 » • • • ,W m llX^’ calculate MpM 2 , . . from 

equation (A21), and tpt 2 , . . . ^tpilX * rom equation (A32). This is programmed in 
FORTRAN subroutine CORODE. 


35 



n n 


APPENDIX B 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


BASIC COREST FORTRAN PROGRAM 


**** PRO GR AH ANALYSIS OF METAL OXIDATION/SCALE VAPORIZATION BY 
MEANS OF PAR AL INEAR KINETIC MODEL • REFERENCE ~ O.WAJZEL, 
J.ELECTROCHEM. SOC* , VOL . I 1 0, NO • 6 , JUNE 1 96 3, PGS .5 04-507 . 

THE PROGRAM PROVIDES THE HISTORY OF THE NET SAMPLE WEIGHT, 
ALLOY METAL CORRODED, AND SCALE THICKNESS. 

**** ADDITIONAL RESULTS PROVIDEO FROM THE BASIC SOLUTION ARE — 

4 1 I CROSSOVER TIME FOR NET SPECIMAN WEIGHT CHANGE OF ZERO 
4 2 1 MAX I MUM NET SAMPLE WEIGHT AND ITS CORRESPONDING TIME; 
4 3 1 WEIGHT OF METAL CORRODED UP TO CROSSOVER TIME; 

19 1EOUILIBRIUM WEI6HT OF SCALE RETAINED ON THE SAMPLE 
AT VERY LARGE TIMES; 

451 ELAPSED TIME FOR AN EXPECTED WEIGHT LOSS OF ALLOY; 

16 1 TIME TO REACH SPECIFIED FRACTIONS OF EQUILIBRIUM 
SCALE THICKNESS; 

471PARALINEAR PARAMETERS ARE DETERMINED FROM THE 

EXPERIMENTAL DATA IN THOSE CASES WHERE THEY ARE 
NOT GIVEN AS INITIAL INPUT. 

**** PROBLEM DEFINITION ANO DATA DESCRIPTION CARDS ARE LISTED AS - 
CARD 1 - FLA6 FOR TYPE OF SCALE LOSS. TRUE rSPALLING, 
FALSE -VAPORIZATION. 

CARD 2 - ALLOY TYPE, ALLOY DENSITY 
CARO 3 - SCALE TYPE, SCALE DENSITY. 

CARD 9 - DATA SOURCE AND RUN DESCRIPTION. 

CARD 5 - EXPERIMENTAL TEMPERATURE ANO PRESSURE. 

CARD 6 - FLAG FOR MICROFILM PLOT REQUEST. TRUE r PLOTS 
WANTED, FALSE = UNWANTED. 

DIMENSION HOURS! lOu 1 , OP 4 1 D 1 , XXX 4101,TP 410 1,TTT <1001 
DIMENSION WW! 1001, TTI 1Q01,XKP4 1001 ,XKVI 1001 
DIMENSION WM4 1 3C1 , WR4 1DD1 , WT4 1 Qul , THKLOS < 1001 ,0X1 THK <1301 
DIMENSION ALLOY 161 , SC ALE 4 6 1 , S OURCE < 8 1 , RUN < 2 1 , TEMP ! 2 1 , PRE S S 4 2 1 

REAL KP,KV 

LOGICAL SPALL, SEMI LG 

COMMON/BLK 1/HOURS, OP, X XX, WW,TT 

COMMON/B LK2/A SMALL , B SMALL , DE NO X ,DENHET 

COMMON / B LK 3/ITYPE,LTVPE,LL,LLX,NUHSET, NUMKP,NUHKV 

COMMON/PAR AMI /ALPHA ,BE TA/PAR AM2/TAUZR0,R A T 10 

COMM ON /C ASE 1 /KP, K V/CASE2/ACA P, BC AP/C A SE 3/WBAR , TB AR , TZERO 

COMMON /€ ASE4/XKP, XK V 

COMMON/ V AL1/WMBAR,WRBAR,WMZER0, WRLIM 

COMMON/S CALE/TP/RUST/TTT 

COMMON/ V OL UP /MM,WR,WT, THKLOS ,OXITHK 

COMMON /INFO/ALLOY, SCALE, SOURCE, RUN, TEMP, PRESS 


• .ARRAYS FOR MICROFILM PLOTTING ROUTINE 

DIMENSION XTITLE 14 1 , YMLEG4 91 , YPLEG 4101,YTLEG491 
COMMON/T ITLES /XTITLE ,YMLEG, YPLEG, YTLEG 
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COMHON/LE6END/XLEGIZOI,YLE6I20I,NX,NY 
COMHON/BLKPLT/KKKI 101 
COMHON/F ILH/SEMILG 


C 

CSM*«M*«**** ****************************************************** 
C *** INITIALIZE DATA AND STOICHIOHE TRIG COEFFICIENTS 

C 

CALL BEGIN 


CS*** ********* ****************************************************** 


C *** DESCRIPTION OF DATA TYPE 

C I TYPE-1 THE DATA SET CONSISTS OF THE CURVE PARA- 

C METER PAIR WB AR, TB AR • 

C ITYPE=2 THE DATA SET CONTAINS NUHSET PAIRED VALUES 

C OF THE PARALINE AR OXIDATION-REACTION PARA- 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


HETER KP, AND THE VAPORIZATION PARAMETER KV* 

I TYPE -3 GIVEN NUHSET PAIRS OF EXPERIMENTAL DATA 

POINTS WU4l>,TT<I»,NUHKP VALUES OF THE PARA- 
LINEAR PARAMETER KP — ILTYPE=1I OR NUMKV 
VALUES OF THE VAP0R12AT ON PARAMETER KV — 
ILTYPE-2 1 • 

I TYPE =9 6IVEN NUHSET PAIRS OF EXPERIMENTAL WEIGHTS 

W1 1 1 1 « W2 II! AND THEIR CORRESPONDING TIMES 
T1<I»,T2<1> . 


CS ********************************************************* ********* 


C *** CYCLE THROUGH DATA VECTOR 

DO 500 1=1, NUHSET 
C 

C *** TEST FOR DATA TYPE 


GO TO 160,70,80,901 ,ITYPE 
C 

60 CON I INUE 


CS************ ****************************************************** 
C **** ITYPErl 

WBARrvWt I » 

TBARrTTII I 

WRITE <6, AO >WBAR,TBAR 
CALL QUERY 1<WBAR,TBAR ) 

CALL OUTPUT 
60 TO 500 
C 
C 

70 CONTINUE 
C 


CS ************ ****************************************************** 
C *** I T YP E=2 

KP=XKP«I I 

KV=XKV<1) 

WRITE 16, *2>KP,KV 
CALL ODER Y 2 I KP , K V ) 

CALL OUTPUT 
60 TO 500 
C 

c 

80 CONTINUE 

CS************ ****************************************************** 
C ♦*** ITYPErT 

C 
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W2Z -U W 1 1 ) 

122 - TTIII 

*** TEST FOR INPUT TYPF OF P ARAL INE AR PARAMETER **** 
60 TO <810,8201, LTYPE 

810 CONTINUE 

DO 811 JKP-1 , NUMKP 
K P — X KP I JK P ) 

WRITE I6,45)KP,WZZ,T2Z 

CALL QUERY 3 I L TY PE , WZZ , IZ7,KP) 

CALL OUTPUT 


811 

CONTINUE 


60 TO SCO 

820 

CONTINUE 


DC 821 JK V ~ 1 , NUMK V 


KV=XKV<JKV) 

WRITE < 6 , *6 )KV,WZZ,TZZ 

CALL QUERY 3ILTYPE,WZZ,IZZ,KV ) 

CALL OUTPUT 
821 CONTINUE 

GO TO 500 

90 CONTINUE 

************* ****************************************************** 
**** ITYPE-M 

N4r2*I 
N 4M IN 1 -N 4-1 
Wlr WWIN4MIN1I 
W2- UW < N4 ) 

Tlr TTIN4HIN1I 
12 - TTIN4 1 
WRITE<6,47IW1,T1 
WRITEI6,48)W2,T2 
CALL QUERY4IW1, W2,T1,T2I 
CALL OUTPUT 
60 TO sro 


500 CONTINUE 

************ ****************************************************** 
**** FORMAT STATEMENTS FOLLOW **** 

4 C FORMATI////10X,26HINPUT DATA WMAX r , F 1 2 . 6 , 1 2 X , 7HTM AX ~ , 

1F12.6,///) 

42 FORMATI ////10X,24HINPUT DATA KP r , E 1 2 . 8 , 1 2 X , 5HK V r , 

1F12.8,///! 

45 FORMATI/ ///1CX,23H1NPUT DATA KP = , F 12 • 8« 1 OX , 7HWV 111 = , 

IF 12.6, 10H MG/Se.CM., 10X f 7HTTI II r f F12.6,6H HOURS///) 

46 FORMAT I ////1QX, 2 3H INPUT DATA KV r , F 1 2 . 8 , 1 0 X , 7HW W < 1) r, 

IF 12.6,1 OH MG/SO. CM. ,10X,7HTTI1) r,F12.6,6H HOURS///) 

47 F ORMA T <////10X,26Hl NPUT DATA WW I 1 ) = , F 1 2 . 6 , 1 5 X , 7HTT <1 I 

1F12.6,6H HOURS//! 

48 FORMAT < 29X,7HWWt2I = , F 1 2 . 6 , 1 $X , 7HT T I 2 ) r t F12.6,6H HOURS//) 

STOP 

END 
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SUBROUTINE BE 6 I N 


C 

C *♦** SUBPROGRAM READS PROGRAM DESCRIPTORS FROM DATA CARDS AND 

C WRIfES INITIAL HE ADIN6S ON OUTPUT. ♦ *** 

C 

LOGICAL SPALL * SE MI LG 

C 

DIMENSION ALLOY 16 1 , SCALE I 6 1 , SOURCE i 8 1 , RUN t 2 1 , T EMP < 2 1 , PRE S S « 7 1 
COMMON /B LK2/ASMALL,BSMALL*0EN0X,0ENMET 
COMMON/B LK 3/1 TYPE t L TYPE, LL , LLX , NUM SE T , NU MK V , N UMK P 
COMMON/P AR AMI /ALPHA t BET A /PAR AM2/TAU2RO, RATIO 
COMMON/ INFO/ ALLOY, SCALE, SOURCE , RUN , T E HP , P R E S S 
COMMON/F ILM/SEMILG 
C 

RE AD I 5 , 1 1SPALL 
READ 15*2 1 ALLOY, DEN MET 
READ 15,3 1SCALE , DENO X 
RE AD 4 5 , A 1S0URCE,RUN 
RE AD t 5 , 5 ITEMP, PRESS 
READ45,6 1SEMILG 
C 

IF 4 SPALL 1 60 TO 20 
WRI TEI6, 1001 
GO TO 30 

20 WRITEI6,2C01 

3 0 WRI TE 16, 3 001 ALLOY, SC ALE 

WRI TE 4 6, GO ITEMP, PRESS 
WRI TE I 6, 5 GDI SOURCE, RUN 
WRI TE46, 6CC)DENMET,DEN0X 
WRI TEI6,73f7l A SMALL, B SMALL 
WRIIE46,80D1ITVPE 
WRI TE46,9GC1 
CALL SfOICH 
CALL 2ER0 

WRI TE46,95D1AL PH A, BETA, RATIO 
C 

1 F0RMATIL51 

2 FORMA? 16 A6,F6.31 

3 FORMAT46 A6,F6.3J 

4 FORMAT ( 8 A6,2A6 1 

5 F0RMAT4 2 A6,2A6 I 

6 F0RMATIL51 
C 

100 FORMAT ( 1H1 , 35X ,62HPARALINEAP KINETICS STUDY OF CORROSION WITH SCAL 
IE VAPORI ZATI0N//1 

200 FORMAT! 1 H1,35X,52HPARALINEAR KINETICS STUDY OF CORROSION WITH SPAL 
1LING//1 

300 FORMAT I 1 HO , 35 X , 11HALL0Y T YPE - , 5X , 6 A6 // 35 X , 1 1 H S C A LE TY PE * , 5 X , 6 A 6 1 
AOO FORMAT! 1 HO , 35 X , 1 A HI EMPER A TURE - , 2 A6 , 1 OX , 1 1 HPRES S U RE = ,2A6» 

500 FORMAT! 1 HO , 35 X , 1 2H0 ATA SOURCE- , 8 A6 , 1 OHRUN NO • r f 2A61 
600 FORMAT! 1HG,20X,18HDENSITIES— ALLOY - ,F6.3,21H GM/CU.CM., SCALE 
1- , F 6 • 3 , 10H GM/CU.CM.l 

700 FORMAT I 1 HO, 20 X ,39HST0I CH10ME TRIG COEFFICIENTS- ASMALL r ,F7.*I,12H 
1 , BSHALL - ,F7.N1 

800 F0RMAT«1H0,2DX,12HDATA TYPE = ,111 

900 FORMAT! 1 BO , 5 X , 1 9HD A T A SPE C IF IC A ? I ON-// 1 5 X , 53HITYPE-1 THE DATA 

I SET CONSISTS OF THE CURVE P AR A-/28X , 22HMETER PAIR WB AR , TBA R .// 

215X,59HI TYPE=2 THE DATA SET CONTAINS NUMKP VALUES OF THE/ 

528 X , A 3 HP ARALINEAR 0 IFF US I ON -RE ACT I ON PARAMETER KP , / 28 X , % 2HA ND NUHK 
*V VALUES OF THE VAPORIZATION P AR A -/28 X , 1 OH METER KV ./ / 1 5 X , 5 2H 1 T YPE 
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5-3 GIVEN NUMSET PAIRS OF EXPERIMENTAL 0 A T A/2 8 X , 4 4HP0I N TS VWIl 

6I,TTIII,NUHKP VALUES OF THE P ARA-/28X , 41 HL INE A R PARAMETER KP — «L 
7T YPE —I I CR NUMKV/2BX , 42HVALUES OF THE VAPORIZATION PARAMETER KV — 

8/28X,llH IL TYPE -21 .// 1 5X , 55HIT VPEr 4 GIVEN NUMSET PAIRS OE EXP 

9ERIMEN T A L WE IGHTS/28X , 4 1HW 1 C 1 1 , W2 1 1 1 AND THEIR CORRESPONDING TIMES 
1/28X,13HT1«U,T2II1 .///I 

950 FORM AT I 1 OX f 8H ALPHA = , FI 0 .9, 12 X , 7HBE ¥ A r , F 10 . 9 , 1 3X , 38HCR0S SOVER T 
1IME/TIME IAT MAX. WEIGHT I r ,F8.6,///1 

C 

RETURN 

END 


SUBROUTINE OUTPUT 
C 

DIMENSION HOURS! I 00 I, QPf 1 G I , XX X 4 1 01 f TP I 1 0 I t TTT 1 100 I 
DIMENSION WWf 10C),TTI 1G01,XKP! luC 1 , X K V I 1 CD 1 

DIMENSION WMI ID01,WR« 1 00 1 , WT 1 1 Qu 1 , THKLOS II CO 1 , OX I THK I 1 QO 1 
DIMENSION ALLOT 16 I * SCALE I 6 l,S0URCE<8l ,RUNI21, TEMP I 21, PRESS I 2) 

C 

REAL KP,KV 

LOGICAL SPALL, SEHILG 
C 

C OHM ON/B LK 1 /HOURS ,GP, XXX, WW, TT 

COMMON/B LK 2/A SMALL, BSMALL, DENOX,DENMET 

COMMON/B LK3/ITYPE,L TYPE, LL,LLX, NUMSET, NUMKP,NUMHV 

COMMON/P ARAM1/ALPHA,BETA/PARAM2/TAUZR0, RATIO 

COMMON/ C ASEI/KP,KV/CASE2/ACAP,BCAP/CASE3/WBAR,TBAR,T2ER0 

COMMON/C ASE4/XKP,XKV 

COMMON/ V ALl/UMBAR,yRB AR , W MZE RO , *RL JH 
COMMON /S CALE/TP/RUST/TTT 
C OMMON /V QLUP/WM,WR,WT, THKL OS , OX ITHK 
COMMON/ I NFO/ ALLOT, SCALE , SOURCE , RUN , TEMP , PR ES S 
C 

C. .NAMED COMMON FOR MICROFILM PLOTTING ROUTINE 

DIMENSION XTITLEI41 , YMLE 6191 ,YPLEG< 101 ,VTLEG|9) 

COMMON/ titles/xtitlf ,ymleg,vpleg,ytleg 

COMMON/B LKPLT/KKK I IQ) 

COMMON/F ILM/SEMILG 
C 

CS»** ********* ********************************************************** 

#** PRINTOUT FOR BULK OF COMPUTED RESULTS *** 

C\ ************ ************* ************* ******************* ****** **** *** 

C 

100 CALL SPEVAL }■ 

WRITE I 6, 12 1ACAP,BCAP,KP,KV, TZERO 1 

WRITE 16, 13IWBAR,WRLIM 

WRITE 16, 14 1TB AR f WMZERO T 

WRITE I 6, 15 JWHBAR,WRBAR 
WRITE16,161 
CALL PROVAL 

WRITE 16, 171 <OP<Kl,TPtKl ,K=1,31 
WR I TE 16,181 
CALL CORODE 
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WRITE46,1914XXX(II!,TT14II!,II=1,LLXI 
WRITE 46*21 I 
CALL CURVE 

WRITE 4 6* 23! 4 HOURS 4 I J! , VM 4 I J ! , WR 4 I J! , WT 4 I J 1 , THKL OS 4 I J ! , 

1 OKI THK 4IJI*IJ=1*LL! 

..REOUEST FOR HICROFILH SEMILOG PLOTS 

IF4SEHIL6ICALL PLOTI04LL! 

12 FORMAT4//10X,4HA = , F 1 0. 7 « 6X , 4HB - , F 10 . 7 , 6X * 5HKP = ,F10.7*5X, 

1SHK V r ,F10.7,5X,17HCR0SS0VER TIME = ,F10.2,6H HOURS///! 

13 FORMAT 4 1 OX , 21HMAX1MUM NET WEIGHT = *F10.5,1?H MG/SO .C H .* 23X , 24HLI 
1MIT1MG SCALE WEI6HT = ,F10.5,11H M6/S0 .CM . /// I 

14 FORMAT (10X.29HTIME AT MAXIMUM NET WEIGHT = ,F10.2,7H HOURS, 14X, 38 
1HALLOY C ORROOEO UP TO CROSSOVER TIME = .F10.4.10H MG/SO. CM. /// 1 

15 FORMAT! 1 OX. 3 7H ALLOY CONSUMED! AT MAX. NET WEIGHT! = .FIO.V.ICH MG/S 
10.CM..3X .35HSCALE WEIGHT I AT MAX.NET WEIGHT! = .F10.4.10H MG/SO. CM. 
2////I 

16 FORMAT 4 4 OX . 20HFRAC TI ON OF LIMITING, 14X, 12HELAPSE0 TIME/44X, 12HSCAL 
IE WEI6HT ,21X, 7H4H0URS ! // ! 

17 FORMATI48X.F6.5.22X.F10.2,/! 

18 FORMAT!/ ///42X, 15H ALLOY CORROSION, 16X , 12HELAPSE0 T IME /43X , 1 1H ( MG/S 
10. CM. I, 2 IX. 7H I HOURS!//! 

19 FORMAT44 5X , F 7 . 2 , 24 X , F 1 0 . 2 , / ! 

21 F ORMAT 4 / /// 1 3 X, 12 HE LAPSED TIME , 5X , 1 4HC0RR0 SI ON LOSS, 6 X , 1 2HS C ALE WE 
1I6HT.5X, 19HNET SPECIMAN WE IGHT , 4X , 1 4HC0R RO S ION LOSS, 4 X , 1 5HSCALE TH 
2ICKNESS / 16X.7H! HOURS | ,9X,11H (MG/SQ.CM. 1, 7X, 11H4MG/ SO. CM. I , 1 CX , 

31 1H! M6/S O.CM.!,6X,13H (MILLIMETERS l,6X, 13H4MILLIMETERS !//! 

2 3 FORMAT! 1 5X ,F10 .2,7X ,F12 .4 » 7 X »F 12. 4 ,9X,F12.4, 14X.F 8.6, 10X ,F8.6, / 1 

RETURN 

END 


SUBROUTINE PL0TI04LINK! 

... PLOT OXIDATION-CORROSION DATA. CURVES OF ALLOY WEIGHT LOSS, 
... SCALE PRODUCT, AND NET WEIGHT LOSS VS. TIME.. 

DIMENSION CHAR I (12! , CH AR 1 ( 6 1 , CHAR 2 ! 8 I 
DIMENSION XTI TLE44 I ,YMLE6!9! , YPLEG ! 10) , YTLEG4 9 I 

DIMENSION ALLOY 161, SCALE 16!, SOURCE 48! .RUN I 21, TEMP! 21, PRESS! 21 

COMNON/B LK 1 /HOURS 4100!, OP (1C!,XXX41C!,WM(100I«TT(100! 

COMMON/ VOLUP/WM HDD l.WRl 100 1 , W T ( 1 OQ 1 , THKLOS 4 1 00 1 , OXI THK 4 100 1 

COMMON/ 1 NFO/ ALLOY, SCALE, SOURCE, RUN, TEMP, PRESS 
COMMON/ TITLES /XTI TLE , YMLEG , YPLEG , YTLE6 
COMMON/LE6END/XLEG4 20 1 , YLEG 4 20 ! ,NX , NY 
C0MM0N/BLKPLT/KMK4 10! 

... LOAD LEGEND ARRAYS FROM COMMON BLOCK INPUT 
DO 10 1 = 1*8 

CHAR 2 41! -SOURCE 4 11 
J=I*6 
K =1 2 
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L — 1*4 

IFII.6T.2I GO TO 
CHARI «II=TEMPII ) 

CHARI IKI-PRESST 1 1 
CHARI (L ) —RUN til 
5 IF 1 1 . GT .6 1 GO TO 10 

CHART ill 3 ALLOY III 
CHART «J)=SCALE(I I 
10 CONTINUE 

C 

C. .OVERWRITE KKM61 BY SUBROUTINE CALLING ARGUMENT 
KKK I 61 -L INK 

C 

C. .ABSCISSA TITLE ARRAY 
DO 20 IX -1 * 4 

XLEGI IX I -X T I TLE 1 1 X I 
20 CONTINUE 

NX-19 

C 

C..SEMILOG PLOT OF ALLOY CORROSION HISTORY 
00 30 IM-1,9 
YLEGI IMjrVMLEGTIMI 
30 CONTINUE 

NY- 5 2 

CALL LRTLEGICHART,72) 

CALL LRLE6NICHAR1. 36.0,1. ,7., C.I 
CALL LRLEGNICHAR2.48, 3, l.,6.5,G. I 
CALL LCGXYIH0URS,WM,2,KKK I 

C 

C . .SEMILOG PLOT OF PRODUCT SCALE HISTORY 
00 40 IP-1,10 
YLEGI IP J = YP LEG TIP! 

40 CONTINUE 

N Y-6 P 

CALL LRTLEGICKART, 721 
CALL LRLEGN I CHAR 1.36,0,1., 7., 0.1 
CALL LRLEGNTCHAR2, 48,0,1., 6. 5«u.l 
CALL L0GXY(H0URS,WR,2,KKKI 

C 

C.. SEMILOG PLOT OF NET SAMPLE WEIGHT HISTORY 
DO 50 IT-1,9 
YLEGI IT»=YTLEGIITI 
50 CONTINUE 

NYrSO 

CALL LR TLEG I CHART, 721 

CALL LRLEGNICHAR1,36,0,1.,7.,G.I 

CALL LRLE6NICHAR2,48,C,1., 6.5,0. » 

CALL LOG XY(H0URS,WT,2,KKKI 

RETURN 

END 


SUBROUTINE STOICH 

COMMON/B LK2/A SMALL , B SMALL , DENO X ,DENHE T 
COMMON/P ARAM1 /ALPHA ,BET A 
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c 

Cl ************ ****************************************************** 

c 

C *** FUNCTIONS OF THE STOICHI OHETRI C COEFFICIENT ASHALL *** 

C ALPHA - DIHENS IONLESS HAXIHUH NET SAHPLE WEIGHT 

C BETA - D I HENS IONLESS TIHE AT NET SAHPLE WEIGHT ALPHA 

C 

c ************* ****************************************************** 

c 

U=8 SHALL/ ASHALL 
XL“ ALOG I U I 
ALPHA =U-1 .-XL 
BE T A- XL -I 1 ♦/B SHALL 1 
RE TURN 
C 

ENO 


SUBROUTINE ZERO 
C 

COHHON/B LK 2/ A SMALL * B SHALL t DENOX*DFNHET 
COMMON /P ARAH1 /ALPHA f BET A /PAR AH 2/TAUZRO, RATIO 
C 

Cl ************ ****************************************************** 

c 

c *** functions of the stoichioheteric coefficient asmall *** 

C TAUZRO - PREDICTED DIMENSIONLESS TIHEIFROM PARALINE AR 

C THEORY I AT CROSSOVER POINT FOR SAHPLE WEIGHT 

C RATIO - CROSSOVER TIHE TO TIHE AT MAXIMUM SAHPLE WEIGHT 

C 

cs ************ ****************************************************** 

C 

c 

1 FORMAT 1 1 HQ * 20 X f 4 1 HERR OR ALERT MESSAGE FROM SUBROUTINE ZERO /1H , 

1 25X $ 6 7H* *THE NUMBER OF ITERATIONS EXCEED IDO AND EXITS WITH THESE 
2VALUES**/! 

2 FORMAT! 1H f 3 OX , I6H T AUZRO* A SM ALL r ,F9.5/1H ,3QX f 37HTHE TRIAL ITERA 
I T ION VALUE ♦♦TEST#* - .F9.5//I 

C 

c 

C ** I N I T I A LI ZE FOR FIRST ITERATION ** 

UrBSMALL/ASMALL 
F IRS T — 1 • /BSMALL 
P%l— SORT I FIRST 1 
C 

C ** ITERATE ON NEXT ESTIMATE ** 

DO 55 1=1,100 

OELT A =1 P J 

P JJ=P»I* I ALOG (DELTA 1*U*P J 1/ 1 1 ./DELTA-UI 

ZFF=«PJJ-PJI/PUJ 

TEST=ABS I ZFF ) 

IFITEST.LT. .000011 60 TO 105 
PJ=PJJ 

55 CONTINUE 

WR I TE 1 6 « I ) 

WRI TE 16*2 ) P J J* TEST 
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c 

105 TAyZRO^P JJ/A SMALL 

RATl Q-TAUZRO/BETA 
RETURN 
C 

END 


SUBROUTINE OUERYl(y,T) 

REAL KP,KV 

COHHON/B LK 3/1 T YPE * L TYPE » LL * LLX * nUMSE T «NUHKP*NUKK V 
COHHON/B LK 2/A SM ALL #BSH ALLUDE NO X f DENHEt 
COHHON/P ARAHl ✓ ALPHA t BET A /P*R AH2/TAUZR0, RATIO 
COHHON/CASEl/NP,KV/CASE2/ACAPtBCAP/CASE3/WBAR^TBAR f TZERO 

£ 

CS ************************************************************ ****** 

c 

C *** SUBROUTINE DETERMINES NUMERICAL VALUES OF THE PARALINEAR PARAMETERS 
C FROM ONE OF SEVERAL POSSIBLE TYPES OF INPUT DATA ** 

C 

C *** INPUT AR 6UHENTS WBAR,TBAR FOR ITVPE - 1 *** 

C 

€!♦♦* ********* ****************************************************** 

C 

40 WB A R-W 
TBAR^T 
C 

BCAPryBAR*BETA/< TBAR*ALPHA| 

AC A P-BCAP*WB AR/ALPHA 
AA- ACAP/BCAP 
NP-2»*ACAP/ASHALL**2 
KVrBCAP*B SHALL/ A SMALL 
TZERO -RATIO* TEAR 
RETURN 
C 

END 


SUBROUTINE 0UERY2IP f VI 
REAL KP v KV 

COMHON/C ASEI/KP,KV/CASE2/ACAP*BCAP/CASE3/H BAR* TBAR f TZERO 
COMMON/P ARAM 1/ ALPHA f BETA/PARAM2/ T AUZRO *R ATIO 
COMHON/B LK 2/A SMALL « BSM ALL f OENOX tDENKE T 
COMHON/B LK3/ 1 T YPE* LTYPE t LLf LLX « NUMSET * NUMKP*NUWK V 
C 

CS ************ ****************************************************** 

c 

C *** SUBROUTINE DETERMINES NUMERICAL VALUES OF THE PARALINEAR PARAMETERS 
C FROM ONE OF SEVERAL POSSIBLE TYPES OF INPUT DATA ** 

C 

C *** INPUT ARGUMENTS KP AND KV FOR ITYPE " 2 *** 

C 


44 



u u 


cs ************ ****************************************************** 

c 

60 KP-P 

K¥-V 

C 

ACAPrKP*ASMALL**2/2, 

bcap>kv*a small /bshall 

AA- ACAp/BCAP 
WR A R^AA^ ALPHA 
TBAR-AA*B^TA/BCAP 
T?ER0rRATI0*TBAR 
RETURN 
C 

END 


SUBROUTINE OUE RV3 1 W « T v X ) 

REAL KP,KV 

COHMON/8LK 2/A SHALL, BSHALL ,DENOX t DENMET 

COMHON/BLK 3/1 TYPE, LTYPE, LL,LLX f NUH SE T , NUHK P, NU *K V 

C OMM ON/P AR AMI /ALPHA, BET A/P ARAM2/TAUZR0, RATIO 

COMMON ✓€ ASE1/KP,KV/CASE2/ACAP,BCAP/CASE3/UBAR, TBAR^TZERO 

c 

CS ****************************************************************** 

c 

C *** SUBROUTINE DETERMINES NUMERICAL VALUES OF THE PAR ALlNEAR PARAMETERS 
C FROM ONE OF SEVERAL POSSIBLE TYPES OF INPUT DATA ** 

C 

C *** INPUT ARGUMENTS ONE WEI6HT-TIME DATA PAIR PLUS EITHER HP 

C OR KV FOR ITVPE = 3 ** 

C 

C%*** ********* ****************************************************** 

c 

1 F ORMAT 1 1 HO 9 20 X , A3HERR0R ALERT MESSAGE FROM SUBROUTINE QUERY 3 /1H , 

1 25X , 67H* ♦THE NUMBER OF ITERATIONS EXCEED 100 AnO EXITS WITH THESE 
2VALUES**/! 

2 F ORMAT I 1 H , 30X , 25HTHE PARAMETER **BCAP** = ,E8,4/lH , 30X , 37HTHE TR 
UAL ITERATION VALUE ♦♦TEST** r ,F 9.5///I 

3 FORMAT I 1 H ,30X,27HTHE ABSCISSA TERM **TH** - ,E 8.9/1H ,30X,37HTHE 
1TR1AL ITERATION VALUE **TEST** - ,F9,5///* 

C 

<1 W U 1 - u 

5 TT1“T 

IFILTYPE .EO. 21 GO TO IQ 
C 

C ** COEFFICIENT KP GIVEN ** 

KP^X 

AC A P-KP* A SMALL** 2/2 • 

C 

A A- ACAP/ A SMALL 

BHA Xr |- WW 14S0RTI UU1**2^A.*TT 1/ASMALL I */2.*TTl 

** FIRST GUESS OF BCAP VALUE ** 

6 BlrBMAX 
C 
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00 8 jri,iao 

T1=TT1*B1^*2/ACAP 

P— RETRNO I T 1 1 — T I 

PP-I 1 P M l 2«*Tl/<P + Bl I 

F rTT 1 *B 1**2+NW1 1 — A A *P 

0F=2.*TT1^B1*WW1-AA#PP 

B 2-B 1 -F/DF 

TES T - ABS • IB2-B1 1/B2I 

IFITEST.LT. .0001 1 60 TO 9 

Bl=82 

8 CONTINUE 
HR I TE«6 f 1 I 

HRI TEI6 f 2IB2, TEST 

9 BCAP-B2 
C 

KV=BS«ALL*BCAP/ASMaLL 
UBAR-ACA p*alpha/bc a p 
T B A R- UR A R *B ET A / 4BCAP* ALPHA ) 

T2E RO-TAUZRO+TBAR/BETA 
RETURN 

♦♦COEFFICIENT KV GIVEN ♦♦ 

10 K V - X 

BCA P=KV* A SHALL /B SHALL 

♦♦ ITERATION USING KODIFIED SECANT HETHOD ♦♦ 
CON ST^BS HALL* A SHALL ♦UHl/l TTl+BCAPI 
TR= l./f CONST-1 .1 
YR^EXPI-CONST^TRl 
TL-1. /CONST 

DO 16 I - 1 * 100 

EL=CONST^TL 
YL^E XP I -EL ) 

♦♦ CALCULATE CHORD SLOPE ♦* 
CS“IYR-YL)/ITR-TLI 

♦♦ ESTIMATE MEAN ABSCISSA VALUE ♦♦ 

TH - I 1 .-YL*CS^TLl/fCONSI*CS-l . I 
TEST=ABSITM-TL1 
IFI TEST.LT. •ODOOU GO TO 17 
TL-TM 

16 CONTINUE 
URITEI6, 1 I 
WRITE I6*31TM,TEST 

C 

17 A C A P- TT1^BCAP^^2/TM 
KP“2.*ACA P/A SM ALL **2 
HB A R- AC A ALPHA/6 CAP 

TB A RrHB A R^BETA/fBCA P^ ALPHA l 
TZE RO-TAUZRO+TBAR/BETA 
RE T URN 
C 

END 
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SUBROUTINE OUE RY4 I W I * W2 y T I , T2 I 
REAL KP * KV 

COMMON/B LK 2/ A SMALL v B SHALL y OENO X* OENMET 
COHHON/BLK3/ITTPEyLTYPEyLLyLLXyNUHSETyNUHKP,NUHKV 
COHMON/P ARAH1/ ALPHA y BETA/PARAH2/ T AU2R0y RATIO 
COHHON/C ASEI/KPyKV/CASE2/ACAPyBCAP/CASE3/WBARy TBARyTZERO 
C 

c* ♦♦*****#**♦* **** ****** ************************************** ****** 

c 

C *** SUBROUTINE DETERMINES NUMERICAL VALUES OF THE PARALlNEAR PARAMETERS 
C FROM ONE OF SEVERAL POSSIBLE TYPES OF INPUT DATA ** 

C 

c *** INPUT ARGUMENTS TWO PAIRS OF EXPERIMENTAL WEIGHT-TIME DATA 

C FOR ITYPE - 4 ** 

C 

cs ******** ************************ * ************************* ******** 

c 

1 FORMAT! 1 HOyZOXy 43HERROR ALERT MESSAGE FROM SUBROUTINE QUERY A /1H y 
125X,67«* *THE NUMBER OF ITERATIONS EXCEED 200 ANO EXITS WITH THESE 
2VALUES**/! 

2 FORMAT I 1 H 9 30X y 2THTHE ABSCISSA TERM **X2** - % E8 • A / 1H y 30Xy 3THTHE 
1TRIAL ITERATION VALUE **TEST#* - yF9.5///l 

3 FORMAT! lH0 f 20X f A3HERR0R ALERT MESSAGE FROM SUBROUTINE 0UERY4 /1H y 
125X y 69H* *DI VI S OR IN ITERATION STEP IS ZERO « ERROR RETURN WITH THES 
2E VALUES **/ 1 

4 FORMATI IHOy 20Xy A3HERR0R ALERT MESSAGE FROM SUBROUTINE 0UERY4 /1H f 
125X y 60H* ♦PARA METER C3 IS TOO SMALL * CHOOSE OIFFERENT WEIGHT PAIRS# 
2*/l 

C 

wwl=wx 

WV2-W2 

c 

WlABSrAB SI Wl 1 
WZABS=ABS! W2I 
C 

C ** TEST FIRST WEIGHT FOR CROSSOVER VALUE *** 

IF! W1ABS.LT*. 0001 lGO TO 1001 
C 

C ** TEST SECOND WEIGHT FOR CROSSOVER VALUE *** 

IF! W2ABS.LT. .0001 1G0 TO 1002 
GO TO 80 
C 

C ** FIRST WEI6HT AS CROSSOVER VALUE ♦* 

1001 TZERO-TI 

TB A RrBETA*TZERO/TAUZRO 
T AU 2- BET A*T2/TBAR 

5 EM 2 -RETRNO I T AU2 I 

6 EN2-!EM2— B SMALL* TAU2 I/ASHALL 
WB A RrALPHA#W2/EN2 

GO TO 112 
C 

C *♦ SECOND WEIGHT AS CROSSOVER VALUE *♦ 

1002 TZE R0-T2 

TBA RrBETA*TZERO/TAUZRO 
TAU 1-BE T A*T1/ TB AR 

7 EMI -RETRNO <TAU 1 1 

8 EN 1 -f E Ml — BSMALL*TAU1 I /ASM ALL 
WBA RrALPHA#Wl/ENl 

60 TO 112 
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C 

C ** ROUTINE for both weight data nonzero ** 

c 

80 IEND-200 iMAX NUMBER OF ITERATIONS 

EPSri.E-7 fiTOLERANCE BETWEEN ITERATES 

IER- Q iERROR CODE , ,N0 ERROR, , 

TT 1 -T 1 
TT 2 -T2 

TR A T-TT2/TT 1 
C1=W1/W2 
C2- 1RAT*C 1 
C3-BSMALL*fC2-l. 1 
C 

IFt ABSIC3I .6T..0D0UG0 TO 90 8TEST ON CHOICE OF DATA 

WRITE(6,4J 3ERR0R MESSAGE 

RETURN 
C 

c 

C. .PREPARE ITERATION ON WEGSTEIN SECANT METHOD • 

C 

90 XSTART“TAUZRO 3INI TI AL GUESS FOR ROOT VALUE 

X2-X STAR T 
X 0— X 2/T R AT 
EMJ- RE T R NO I X 0 I 
EH2-RETRN0IX21 

X“IC 1 *EM 2 -EM 1 >/C3 aiMPLIClT RELATION XrF|X) 

A-X-XO 
B r ~A 

XO-X 

X 2-X 0*TR AT 

EMlrR£TRNO«XO> 

EM2-RE TR NO t X2 1 
TES T “X — • C1*EMZ— EH1I/C3 
C 

c 

C. .START ITERATION LOOP 
DO 106 I’ltlEND 
IF < TES T I 101,111,101 
C 

C.. EQUATION IS NOT SATISFIED BY X 

101 B-B/TEST-1. 

IFIB 1102 ,114, 102 

C 

C. .ITERATION IS POSSIBLE 

102 A=a/B iSECANT CORRECTION TERM 

X - X ♦ A 
B-TeST 

xo-x 

X2-X0^TR AT 
EMI -RE TR NO f XQ I 
EM2-RE TR NO I X2 I 
TEST -X— 1 C1*EM2-EM1 J/C3 
C 

C..TEST ON SATISFACTORY ACCURACY 
TOLrEPS 
DrABSIX I 

IFID-1. I 104,104,103 
1U3 TOL“TOL*D 

104 IFIABSIA I-T0LI1C5, 105,106 8FIRST TOLERANCE TEST 

105 IF«ABS«TEST»^1C.*T0L1111,111,1Q6 iSECOND TOLERANCE TEST 
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106 C0NV1NUE 
..END OF ITERATION LOOP 

WRITE (6,1) 

WR I TE 16, ?)X2* TEST 

IER -1 8ERR0R CODE ,,N0 CONVERGENCE AFTER 

RETURN 8IEND ITERATIONS. ERROR RETURN,, 

** CALCULATE UBAR AND TBAR FIRST ** 

111 TBA RrBETA*TT2/X2 

EN2-IEM2-BSMALL*X2I/ASNALL 
WB A RrALPHA*WW2/EN2 

112 BCA PrWB AR*BETA/( TBAR* ALPHA ) 

ACA P=WBAR*BCAP/ALPHA 

KP=2.*ACAP/ASMALL**2 

KV=BSMALL*BCAP/ASNALL 

TZERO-TAUZRO*TBAR/BETA 

RETURN 

1 1 A IER=2 8ERR0R CODE ,, DIVISOR VALUE ZERO,, 

WR I TE ( 6 , 3 ) 

WRITEI6, 21X2, TEST 
RETURN 


END 


SUBROUTINE SPEVAL 

COMMON /B LK2/A SHALL , BSMALL, DENOX , DENHET /P AR AM2/T AUZ RO , RATIO 
COHN ON/ CASE 2/ ACAP ,BC A P/ V AL1/WMBAR « WRB AR , WMZERO ,WRL IM 

************* **♦♦♦♦*♦♦ ************************************ ********* 

*** CALCULATES PHYSICAL WEIGHTS AT SPECIAL TIMES OF INTEREST — 

1. ALLOY CORRODED AT MAXIMUM NET SPECIMAN WEIGHT, AND WHEN NET 

GAIN IS ZERO! THAT IS, THE ZERO CROSSOVER TIME I 

2. LIMITING SCALE WEIGHT AT INFINITE TIME, AND AT MAXIMUM NET 

VE16HT. 

5 ** ********** ****************************************************** 


** INITIALIZE RATIOS OF SYSTEM PARAMETERS ** 

U= BSMALL/ASMALL 
AA= ACAP/BCAP 

♦ ♦RETURN DIME N SI ONAL (PHYS ICAL ) WEIGHTS AT SPECIAL CONDITIONS 
WMZERO- B SMALL* A A*T AUZRO 
WMB AR- AA*ALOGIU) 

WRB AR— AA/ASMALL 
WRL IM — U* A A 
RETURN 
C 

END 
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SUBROUTINE CORODE 

DIMENSION HOURS! 1QC1, OP I 1 01, TP 4 1GI , XXX 1 1 00 I • TTTi 1001 
DIMENSION WWI1Q0I,TT4100I 

COMMON/BEK1 /HOURS* OR, XXX, WW*TT/SC ALE/TP 
COMM0N/B in 3/1 TYPE, L TYPE, LL,LLX,NUM SET, NUHKP,NUMKV 
COMMON/R ARAM I /ALPHA ,BETA /CASE3/W BAR, TBAR,TZERO /RUST /TTT 
C 

C5 ************ ************************************ c ********* ********* 

C 

C *** SUBROUTINE COMPUTES ELAPSED TIME TTT WHEN XXX MILLIGRAMS PER 
C SQUARE CENTIMETER OF ALLOY HAVE BEEN CONSUMED DUE TO 

C CORROSION *** 

C 

Cl ************ ****************************************************** 


c 

c 

c 

C 

c 

c 

c 

c 


5 


c 


** BE6IN LOOP THROUGH INPUT SEQUENCE OF WEIGHT LOSSES ** 

DO 5 I- 1 , LLX 

** TRANSFORM REAL WEIGHT TO DIMENSIONLESS VARIABLE 

WHXXX,THEN DETERMINE CORRESPONDING DIMENSIONLESS 
TIME TAUXXX, AND TRANSFORM IT BACK TO REAL 
TIME TTT *♦ 

WHXXX=ALPHA*XXX4I l/WBAR 
TAUXXX- RE TRNTIWHXXX) 

TTTI II -TAUXX X*TB AR/BE TA 

CONTINUE 

RETURN 


END 


SUBROUTINE CURVE 

DIMENSION H0URS<10CI*QP< 1 01 , TP 4 1 0 1 , XX X I 1 00 ) , TT T 4 2 00 I , WH f 1 CO I 
DIMENSION WW 4 1 00! , TT I 100 I , WT I 1 00 » , WR| 1 OC 1 , THKLOS I 1 0 D I , OX I THM 1 0 C I 
COMMON /BLK i /HOURS , GP , X X X , W W , TT /S C ALE /TP/ RUST/T TT 
C0MM0N/BLK2/A SMALL , BSMALL ,DENOX , DENME T 
COMMON/B LK 3 /I Ty PE ,L TYPE, LL, LLX, NUMSE T , NUMK P , NUMK V 
COMMON/YOLUP/WM,WR,WT,THKLOS,OXITHK/CASE2/ACAP t BCAP 

c 

cs ************ ****************************************************** 

c 

C *** THE PRINCIPAL SUBPROGRAM WHICH COMPUTES THE NET-WEIGHT, TIME 
C PAIRS OE DATA WHICH APPROXIMATES THE EXPERIMENTAL 

C CURVE ACCORDING TO THE PARALINEAR THEORY. AT EACH 

C SPECIFIED TIME THE PROGRAM CALCULATES 

C 1. NET SPECIMAN WEIGHT 

C 2. METAL WEIGHT CONSUMED 

C 3. WEIGHT OF SCALE FORMED 

C 4. LOSS IN METAL THICKNESS 

C 5. THICKNESS OF SCALE LAyER 

C 

Ci ********************************************************* ********* 

C 

c 
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c 


** INITIALIZE THE PARALINEAR COEFFICIENTS ** 
AA^A CAP/BCAP 


C 

C 


c 

c 

c 


c 

c 

c 


c 

5Q0 


c 


END 


** SOLVE DIMENSIONLESS OM-LQUATION AT LL TIMES 
00 SOD J- 1 *LL 

TAU-HOURS 4 J1*BCAP/AA 

OMZ retrnoataui 

** THREE VECTORS STORE PHYSICAL WEIGHTS OF ALLOY CORRODED, 
SCALE FORMED* AND NET SPECIHAN WEIGHT ** 

WMt J I “OH* A A 

WR t J I r A A*B SM ALL* ! QM-TAUJ/ASMALL 
WTt JJ-WR! JI-WM! Jl 

** TWO VECTORS STORE PHYSICAL DEPTHS OF ALLOY CORROSION 
AND SCALE FORMED ** 

THKLOS !JJ— *01*WM IJ1/DENMET 
OXITHK! J I -.01* UR I Jl/DENOX 

continue 

RETURN 


FUNCTION RETRNOITAUI 
C 

cs*** ********* ********** ************************************** ****** 

C 

C **♦ SUBPROGRAM RETURNS VALUE OF DIMENSIONLESS ALLOY-CORROSION VARIABLE 
C OM AT TIME TAU ** 

C 

CS** ************** ************* **************************** ********* 

c 

1 FORMAT! 1H0,20X«43HERR0R ALERT MESSAGE FROM SUBROUTINE RETRNO /1H , 
125X,67h*YTHE NUMBER OF ITERATIONS EXCEED 1 00 AND EXITS WITH THESE 
2VALUES**/! 

2 FORMAT! 1H ,30X,3«»HTHE ALLOY CORROSION TERM **0M** 3 ,F9.5/1H ,30X f 
137HTHE TRIAL ITERATION VALUE **TEST** - ,F9.5//1 

C 

C • # TEST FOR SMALL VALUES OF DIMENSIONLESS TIME.# 

XrTAU 

2FIX.6T. .01 160 TO ID 
C 

C . .THIS SECTION USES THE NE WTON—R A PH SON METHOD FOR ROOT SOLUTION OF A 
C.. TRUNCATED POWER SERIES APPROXIMATION OF THE EXPONENTIAL AT VERY SMALL 
C..TIMES... 

C 

C. .FIRST GUESS AT ROOT POSITION 
PHI— 2 . 

DO 3 LLR ”1,100 
PHI X-PHI *X 

FUNC=1.-I . 5-SORT! PHIXI/6. ♦PHI X/2*. I *PHI 
C 

C. .NEWTON-RAPHSON 

DFUNC-— • 5 ♦ .25* SORT I PH IX 1 — PHIX/12 • 

DEL PHI-— FUnC/DFUNC 

TES TrABS! DELPHI/PHI J 

IFI TEST.LT..QDODD1BGO TO 5 
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PHI-PHI+DELPHI 
3 CONTINUE 
URITEI6, II 
Y 2— S OR T < PH IX | 

WRITE! 6 « 2IY2, TEST 
5 RETRNO-S QRT4 PHI X I 
RETURN 

** THIS SECTION USES NEW TON— R APHSON TO SOLVE COHPLETE PAR AL INEAR 
EQUATION TOR DIKENSINLESS CORROSION WEIGHT .. 


*♦ SET UPPER BOUND ON THE DEPENDENT VARIABLE 
10 C=1.*TAU 


C 


c 


** INITIALIZE FOR FIRST ITERATION ** 
Ylr C 


11 

12 

13 


*♦ BEGIN NEWTONS METHOD** 

DO 13 J-1,1Q0 

PHIr Y 1 ♦EXP! — Y 1 ) — C 
TEST = ABSIPHI1 
IF ITEST.LT ••000011 GO TO 1* 
OPHIr 1 *-E XP! -Y 1 1 
Y2 - Y 1 - PHI/DPHI 
Ylr Y2 
CONTINUE 
WRITE 16 t II 
WR 1 TE f 6 1 2 1 Y2 » TE ST 


1A RE T RNQ- Y 1 

RETURN 


END 


FUNCTION RETRNTfXI 
C 

Cs ************ ****************************************************** 

C 

C *** SUBPROGRAM EXPLICITLY COMPUTES DIMENSIONLESS TIME FROM THE 
C OM EOUATION *** 

C 

cs*** ********* ****************************************************** 

C 

RE TRNT-! X-l . I ♦£ XP I -X I 

RETURN 

END 
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APPENDIX C 


ACAP 

ALPHA 

ASMALL 

BCAP 

BETA 

BSMALL 

DENMET 

DENOX 

KP 

KV 

NUMKP 

NUMKV 

NUMSET 

OXITHK 

RATIO 

TAU 

TAUZRO 

TBAR 

THKLOS 

WBAR 

WRLIM 

WRMAX 

XKP 

XKV 


KEY FORTRAN SYMBOLS FOR COREST 

parabolic scaling constant, A, in terms of W m 
dimensionless maximum net specific sample weight 
stoichiometric constant, a 

linear-vaporization constant, B, in terms of W m 

dimensionless time at net specific sample weight ALPHA 

stoichiometric constant, b 

density of the metal consumed, P m 

density of the retained oxide, P r 

parabolic scaling constant, kp 

linear vaporization constant, 

number of kp- values in input array XKP 

number of ky- values in input array XKV 

number of cases in a COREST run 

computed oxide thickness, mm 

ratio of crossover time to T, tg/t 

dimensionless time, B t/A 

computed dimensionless crossover time 

time to reach maximum specific sample weight change, t, hr 

computed metal loss, mm 

maximum specific sample weight change, AW 

limiting retained specific scale weight 

corresponding retained specific scale weight at AW 

input array of kp- values for data types 1, 2, and 3 

input array of k y - values for data types 1, 2, and 3 
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APPENDIX D 


TYPICAL OUTPUT LISTING FOR COREST 


BLCCK DATA 

..DATA AND ARRAyS FOR PL0TT1N6 SUBROUTINE PLOTIOILLJ 
01 PENS ION XTITLE(4I » YMLEG 19 I * YPLE6 1101 , YTLEG 191 
CO PMON/TITLES/X TITLE, VMLE6,YPLE6, YTLEG 
COPMON/BLKPLT/KKKI 101 
DA1A KKK/3,-1 ,1, 0,0,1, **0/ 

0 A 1A XTI TLE/1 9HEXP0SURE TIME, HOURS/ 

DATA YMLEG/52HWEI6HT OF CORRODED ALLOY, ML6MS PER SQUARE CENTIMETER 
1/ 

OATA YP LEG/60 HWEI6HT OF ADHERING PRODUCT SCALE, MLGMS PER SQUARE CE 
INT IMETER/ 

data ytleg/sohnet sample weight gain, mlgms per square centimeter/ 


..ARRAYS AND DATA FOR PROBLEM SOLUTION 

01 PENSION HOURS 1 1001,GPI 101 , XXX ( 10 I, ITT (100 1 
01 PENSION UW ( 100 1, TT ( 100 1 ,XKP( 100 1,XKV( 100 1 

C0PM0N/BLK1/H0URS,QP,XXX ,*w ,tt/rust/ttt 
COPMON/BLK2/ASMALL,BSMALL,DENOX,DENMET 
C0PH0N/BLK3/ITYPE,LTYPE,LL«LLX,nUH5ET,NUMKV ,NUMKP 
C0PM0N/CASE4/XKP,XKV 
INPUT FOLLOWS 

DATA AS M ALL, BSMALL/2 . 1667,3.1667/ 

OATA IT YpE ,LL ,LLX/4 ,17,3/ 

OATA NUMSET/2/ 

DATA OP/. 5,. 9, .95, 7*0.0/ 

DATA XXX/71.9, 359.5, 719.0, 7*0,0/ 

DATA HOURS/ .31, .69, .86, 1.11 ,1.64,3.06,5.83,20.63,23.61 ,30 .56 , 
* 55. 56, 86. 11, 116. 67,1 66. 67, 194. 99,222. 22, 236. 11, 83*100./ 

DATA WU/2. 50, 3. 29, 3. 24, 3. 13, 96*0.0/ 

DATA IT/ 30. 56,66.11 ,221 .22 ,236 . 1 1 ,96*100 . / 

C 

ENC 
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PARAL1NEAR KINETICS STUDY OF CORROSION UITH SCALE VAPORIZATION 


ALLOY TYPE- CHROMIUM OXIDE FORMERS 2 SETS 

SCALE TYPE- CR203 TO CR03 

TEMPERATURE = 1200 C PRESSURE = 1 ATM. 

DATA SOURCE- CR-.2Y TEDMON IREF.4J BEFORE t AFTER UBAR RUN NO. = 

DENSITIES- ALLOY = 7.190 6M/CU. CM. v SCALE = 5.210 GM/CU.CH. 

STOICHIOMETRIC COEFFICIENTS- ASMALL = 2.1667 v BSMALL = 3.1667 

DATA TYPE = 4 
OATA SPECIFIC ATI ON- 

ITYPE=1 THE DATA SET CONSISTS OF THE CURVE PARA- 

METER PAIR UBARtTBAR . 

ITYPE-2 THE DATA SET CONTAINS NUHKP VALUES OF THE 

PAR ALINEAR DIFFUSION-REACTION PARAMETER KP f 
AND NUHKV VALUES OF THE VAPORIZATION PARA- 
METER KV • 

ITYPE=3 6IVEN NUMSET PAIRS OF EXPERIMENTAL OATA 

POINTS UUlIIftTTtll » NUHK P VALUES OF THE PARA- 
LINEAR PARAMETER KP — ILTYPE=11 OR NUHKV 
VALUES OF THE VAPORIZATION PARAMETER KV -- 
«LTYPE=21 . 

I?YPE=4 6IVEN NUMSET PAIRS OF EXPERIMENTAL WEIGHTS 

U 1 1 1 1 V U2 1 1 1 AND THEIR CORRESPONDING TIMES 
Tlll) v T2m . 


ALPHA = .0820465 94 BETA = .063698616 CROSSOVER T IHE/T1HE 1 AT MAX. WE16HTI = 4.034745 


INPUT DATA UU11! = 2.500000 TT(11 = 30.560000 HOURS 


NW 1 2 1 - 3.290000 TTI2 ) - 66.110000 HOURS 


A = .8367918 B = .0201046 KP = .3564921 KV - .0293835 CROSSOVER TIME = 532.07 HOURS 

CJI 

CJ1 


Ol 

05 


MAX 1 MUM NCI WEIGH 


3.«in93 MG/SQ.CH 


LI HI TINS SCALE WEIGHT 


60.83168 MG/SQ.CM./ 


TIME A I MAXIMUM NET WEIGHT = 131.87 HOURS ALLOY CORRODED UP TO CROSSOVER TIME = 33.8786 MG/SQ.CH. 

ALLOY CONSUMED t A 1 MAX. NET WEIGHT I r 15.7989 MG/SQ.CM. SCALE WEIGHT (AT HAX.NET WEIGHTS = , 19.2098 HG7SQ.CH. 


FRACTION OF LIMITING 
SCALE WEIGHT 


ELAPSED TIME 
(HOURS ) 


.50000 

•90000 

.95000 


399.87 

2903.72 

•1235.21 


ALLOY CORROSION 
(HG/SQ.CM. I 


ELAPSED TIME 
(HOURS ) 


71.90 

359.50 

719.00 


1873.99 

15811.57 

33692.67 


ELAPSED TIME 
1 HOURS 1 

CORROSION LOSS 
IMG/SQ.CM. 1 

SCALE WEIGHT 
I MG/SO . CM « I 

NET SPECIHAN WEIGHT 
IMG/SO. CM.) 

CORROSION LOSS 
1 HILL JHE TEAS) 

SCALE THICKNESS 
•MILLIMETERS) 

.31 

.7228 

1.0867 

.3283 

.001005 

•002009 

.69 

1.0792 

1.5571 

.8778 

•001501 

.002989 

• 86 

1.2055 

1.7366 

.5311 

•001677 

•003333 

1*11 

1.3708 

1.9703 

.5999 

.001906 

•003782 

1.68 

1.6678 

2.3893 

.7215 

.002320 

•008586 

3.06 

2.2837 

3.2878 

• 968) 

•D03176 

•006238 

5*83 

3.1632 

8.85)8 

1*2886 

•008399 

• 008585 

20.83 

6.0873 

8.2263 

2.1790 

.008811 

.0)5789 



23.61 


6.9983 


8.7306 


30.56 

7.3629 

9.8629 

55,56 

10.0308 

13.0278 

86.11 

12.6107 

15.9007 

116.67 

19.8011 

18.2090 

166.67 

17.8971 

21.2598 

199.99 

19.9919 

22.7016 

222.22 

20.8963 

29.0110 

236.11 

21.5999 

29.6238 


INPUT 

DATA 

UUU1 = 

3.290000 

TTU 



UU 1 2 1 = 

3.130000 

TT (2 

A = 

.9907335 

B = 

.0229299 

KP - .9220797 

HAXIHUH NET UEIGRT 

- 3. 

62991 HG/SQ.CH 

I. 

TIHE 

AT HAXIHUH NET HEIGHT 

□ 125.50 

HOURS 

ALLOV 

CONSUMED! A 1 

MAX. NET 

UEIGHTI = 

16.7661 HG/SQ.CH. 


FRACTION OF LIMITING 
SCALE WEIGHT 

.50000 

.90000 

•95000 

CJ1 

<3 


2.2823 

•006968 

.016757 

2.5000 

.010290 

.018930 

2.9970 

•01395] 

.02 5005 

3.2900 

.017539 

.030519 

3. 9030 

•020586 

.039991 

3.3627 

•029892 

.090806 

3.2597 

•027090 

.093573 

3.1197 

.029063 

.096086 

3.0290 

•030035 

.097263 

221.219999 HOURS 



236.110001 HOURS 




KU = .0327739 CROSSOVER TIHE = 506.36 HOURS 

LI MITING SCALE WEIGHT = 69.57212 HG/SQ.CH./ 

ALLOT CORRODED UP TO CROSSOVER TIHE □ 35.9575 HG/SQ.CH. 

SCALE UEI6HTMT HAX.NET UEIGHTI □ , 20.3910 HG/SQ.CH. 


ELAPSED TIHE 
(HOURS ) 

380.59 

2763.92 


9030.56 


cn 

CO 


ALLOY CORROSION 
4HG/SQ.CM.J 

71.90 

359.50 

719*00 


ELAPSED TIME 
♦ HOURS 1 

CORROSION LOSS 
CK6/5G.CM. 1 

SCALE WEIGHT 
♦MC/SO.CH.l 

.31 

.7861 

1.1387 

.69 

1.1745 

3.6939 

• 86 

1.3119 

1*8891 

1.11 

1.4914 

2.1433 

1.64 

1.8150 

2.5989 

3.06 

2.4855 

3.5323 

5.83 

3.4429 

4.8409 

20.83 

6.5841 

8.9401 

23.63 

7.0209 

9.4875 

30 . 56 

8.0169 

10.7154 

55*56 

10.9264 

14.1469 

86.11 

13.7390 

17.2579 

116.67 

16.1290 

19.7493 

166.67 

19.5088 

23.0503 

194.44 

21.1960 

24.6061 

222.22 

22.7849 

26.0178 


236.11 


23.5482 


26.6782 


ELAPSED TIME 
♦ HOURS 1 

1623.14 
14062. 02 
30093.12 


NET SPECTMAN WEIGHT 
♦HG/SO.CH. 1 

CORROSION LOSS 
i MILLIMETERS 1 

SCALE THICKNESS 
♦ MILLIMETERS J 

.3526 

•001093 

.002186 

.5194 

•001633 

.003251 

.5773 

•001825 

.00 3626 

.6519 

.002074 

.004114 

.7839 

.002524 

.004988 

1.0468 

.003457 

•006780 

1.3980 

.004789 

.009292 

2.3561 

.009157 

•017160 

2.4668 

.009765 

•018210 

2.6985 

.011150 

• 020567 

3.2215 

.015195 

.027153 

3.5188 

•019109 

•033124 

3.6203 

.022433 

.037906 

3.5415 

.027133 

.044242 

3.4101 

.029480 

.047229 

3.2329 

.031690 

.049938 

3.1300 

.032751 

.051206 
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TABLE I. - STOICHIOMETRIC FACTOR* 1 

[a - (M-fX 4- M 2 y)/Oz; b = a + 1.] 

(a ) One- meta I oxides (b ) Two- metal oxides 
(l.e M spinels) 


Oxide 

a 

Th° 2 

7.2512 

CoO 

3.6831 

NiO 

3.6694 

FeO 

3.4906 

MnO 

3.4338 

Fe S°4 

2.6180 

FegOg 

2.3271 

Cr 2°3 

2. 1667 

MOOg 

1. 9988 

AlgOg 

L 1242 

t -' r< - ) 3 

! 1.0833 

Si0 2 

.8778 


Oxide 

a 

CoCr 2 04 

2.5458 

NiCr 2 G 4 

2. 5423 

FeCr 9 0 4 

2. 4977 

MnCr 9 0 4 

2.4334 

AtCr.dh 4 

2.0466 

CoAIg0 4 

1. 7639 

NiAlg0 4 

1.7605 

FeAI 2 G 4 

1.7158 

MnAlgC 4 

1.7016 

CrAlg0 4 

1.6556 


'"For example, for Cr 0 CL, a - iMi 00 x 2) + 0 _ 
" 6 16X3 

and for NiCr 9 0 4 , a - 11 * (5-2. OQx 2) 


2, 1667 
- = 2. 5423 


TABLE H. - INPUT DATA FOR PREDOMINANTLY OgOg-FORMING ALLOYS FOR 
PARA LINEAR OXIDATION ANALYSIS (COREST) 

[Scale, CrgO^; stoichiometric constants: a - 2. 1667, b - 3. 1667; density of Or, p m = 7. 19 g/cm 
density of Cr„Oo, p„ - 5.21 g/cm 3 J 


Test conditions 


Before AW 


After AW 



Time, 

hr 

Temper- 

ature, 

°C 

Number 
of data 
points 

Alloy A a 
(ref. 6) 

1280 

1200 

31 

Cr-0.2Y 
(ref. 4) 

240 

1200 

17 

Ni-40Cr 

1020 

1200 

87 

Tophet 30 b 
2-37 

2085 

1200 

55 

Tophet 30° 
2-35 

2688 

1150 

71 

Tophet 30 b 
2-42 

1850 

1150 

41 


Ni-2QCr-3Fe-0, 4Si. 
Ni-29.6Cr-l.43SL 


mum AW in paraHnear oxidation , AW, mg/cm 


t t - 100, AW X =: 3. 21 | t = 2i94 i ± « 780, AW X m. 2. 09 

to - 200, AW 2 -3.59 AW =• 3, 62 tg - 1280, AWg - -0.20 


t x ~ 30.56, AWj * 2.50 t - 166.79 t x ^222.22, AW 
t 2 * 86.11, AWg » 3.29 AW =,3.48 tg -236.11, AW 


t % ^ 100, AW r x - 2. 859 t = 250 t x * 400, AWj = 2. 920 

t~ =• 200, AWg =• 3.204 AW = 3. 240 tg - 900, AWg - 0. 669 


t x = 100, AW X = 1, 955 t =■ 225 


t x » 1200, AW X - 0. 875 


t 2 * 196, AWg =2.249 AW -2.297 tg =2050, AWg = -1. 061 


t x s 40, AW X - 1.670 
tg - 80, AW 2 - 1. 871 



t » 110 t x * 700, AW X 
AW - 1. 957 tg -2300, AW 


t * 420 t x = 990, AW X 
AW f 2. 382 t 9 E 1800, AW, 





















































TABLE m. - DERIVED PARALINE AR ESTIMATES FROM COREST OUTPUT FOR SIX 


Cr 2 Og- FORMING ALLOYS FOR THREE TIME STAGES OF INPUT 


(a) Before AW 


Alloy 

Fit, a 

Parabolic 

Linear 

Limiting value of 
specific weight of 
metal oxide 

Time to reach- 

* O 

img/cm 

scaling 

constant, 

vaporization 

constant, 

*v 

2 

AW, mg/cm 

W m , mg/cm 2 
(Mm) 


formed and re- 




0 

71.9 (100) 

tained, in para- 





linear oxidation, 
W r,lim 

Time, 

t, hr 

Alloy A 

1.191 

0. 2344 

0.01833 

64. 13 

899.11 

2913.91 

(ref. 6) 







Cr-0. 2Y 

.156 

.3565 

.02938 

60. 83 

532.09 

1874. 00 

(ref. 4) 
Ni-40Cr 

.345 

.1850 

. 01620 

57.26 

908. 56 

3516.48 

Tophet 30 

1.841 

. 08013 

. 00984 

40.85 

1067. 46 

6848.28 

2-37 







Tophet 30 

15.635 

.1579 

.02368 

33.44 

362.90 

3086.67 

2-35 







Tophet 30 

3.682 

. 1028 

.01322 

38.99 

758. 16 

5199.20 

2-42 








(b) At AW 


Alloy A 
(ref. 6) 

0.297 

0. 1797 

0. 01397 

64.48 

1186.21 

3 810.58 

Cr-0. 2Y 
(ref. 4) 

.206 

.2929 

. 02369 

61.99 

672.59 

2 299.64 

Ni-40Cr 

.243 

.1693 

. 01471 

57.72 

1008. 69 

3 856. 57 

Tophet 30 
2-37 

2.639 

. 09453 

.01158 

40. 92 

907. 82 

5 810. 11 

Tophet 30 
2-35 

12.426 

. 1404 

. 02019 

34. 86 

443. 82 

3 563.35 

Tophet 30 

.621 

. 05446 

. 00644 

42.43 

1694. 59 

10 289. 80 


(c) After AW 


Alloy A 

0. 469 

0. 1339 

0. 01179 

56.95 

124. 11 

4 845.06 

(ref. 6) 
Cr-0. 2 Y 

. 184 

.4434 

. 03397 

65.45 

495. 17 

1 553.64 

(ref. 4) 
Ni-40Cr 

.270 

. 1578 

.01418 

55.83 

1012. 19 

4 074. 22 

Tophet 30 

.352 

. 03507 

. 00529 

33.21 

1612.24 

13 839. 10 

2-37 







Tophet 30 

.460 

.01935 

.00313 

31.03 

2549. 98 

24 018. 00 

2-35 







Tophet 30 

.360 

. 02698 

. 00359 

37.68 

2697.92 

19 421.54 

2-42 








a Estimated error is averaged over whole curve. 




TABLE IV. - DERIVED PARALINE AR VALUES FOR PREDOMINANTLY CrgOg- FORMING ALLOYS WITH GIVEN PARABOLIC SCALING 


AND LINEAR VAPORIZATION CONSTANTS COMPARED WITH PARABOLIC OXIDATION (k y = 0) 


Tempera- 

ture, 

°C 

Parabolic 
scaling con- 
stant, 

V 

(mg/cm 2 ) 2 /hr 

Linear 

vaporization 
. a 

Metal consumed, W m , mg/cm 2 
(Thickness change, /im) 

Limiting 
specific 
weight 
of Cr 2 Og 
on alloy 
surface, 

W r,lim’ 

mg/cm 2 

Specific sample 
weight change 
(maximum pos- 
itive value), 
AW, 

mg/cm 2 

Time to reach 
maximum AW 
in paralinear 
oxidation AW, 
"t, 
hr 

Time to cross 
zero axis 
'where AW = 0 
mg/cm 2 ), 

t Q , 

hr 

Specific weight 
of Cr consumed 
W m at cross- 
over time tQ, 

mg/cm 2 | 

constant, 

v 

71.9 (100) 

359.5 (500) 

719 (1000) 

2 

mg/cm /hr 

Time to reach specific weight of 
Cr consumed per unit area due 
to oxide buildup and scale 
vaporization W , 
t, hr 

1000 

Low, 0. 002 

0. 00098 

96 788 

525 702 

1 061 845 

10.23 

0.57 

665 

2 684 

5.70 



.0098 

10 618 

53 510 

107 124 

1.02 

.06 

6.65 

26.84 

.57 



0 

550 592 

13 764 810 

55 059 241 







High, 0. 05 

.00098 

19 294 

308 580 

815 542 

255. 82 

14.36 

16 628 

67 088 

142.45 



.0098 

8 155 

51 004 

104 618 

25.58 

1.44 

166.28 

670.88 

14.25 



0 

22 024 

550 589 

2 202 370 






1050 

Low, 0. 0053 

0.0021 

44 016 

244 174 

494 375 

12.65 

0.71 

384 

1 549 

7.05 



.021 

4 944 

24 960 

49 980 

1.27 

.07 

3.84 

15.49 

.70 ! 



0 

207 771 

5 194 268 

20 777 072 







High, 0.130 

.0021 

7 590 

129 593 

357 600 

310.39 

17.42 

9 415 

37 987 

172.84 



.021 

3 576 

23 542 

48 562 

31.04 

1.74 

94.15 

379.87 

17.28 



0 

8 471 

211 766 

847 065 






1100 

Low, 0.0122 

0. 0042 

21 555 

121 632 

246 733 

14.56 

0.82 

221 

, 891 

1 8.11 



.042 

2 467 

12 475 

24 985 

1.46 

.08 

2.21 

8.91 

| .81 



0 

90 261 

2 256 526 

9 026 105 





: 


High, 0.300 

.0042 

3 336 

59 492 

169 463 

358.14 

20.11 

5 432 

21 916 

199.43 



.042 

1 695 

11 657 

24 167 

35.81 

2.01 

54.32 

219.16 

; 19.94 ’ 



0 

3 671 

91 765 

367 061 




i 

i 

1150 

Low, 0. 029 

0. 0082 

10 658 

61 913 

125 989 

17.73 

1.00 

138 

559 

1 9.87 



.082 

1 260 

6 386 

12 794 

1.77 

.10 

1.38 

5.59 

.99 



0 

37 972 

949 267 

3 797 189 





1 

t 


High, 0.700 

.0082 

1 452 

27 172 

80 435 

428.02 

24.03 

3 325 

13 415 

| 238.35 



.082 

804 

5 886 

12 293 

42.80 

2.40 

33.25 

134. 15 

23.83 



0 

1 573 

39 328 

157 312 




1 

j 

1200 

Low, 0. 056 

0. 0145 

5 918 

34 900 

71 136 

19.36 

1.09 

85 

343 

10.78 



.145 

711 

3 610 

7 234 

1.94 

.11 

.85 

3.43 

1.08 



0 

19 664 

491 600 

1 966 401 







High, 1.35 

.0145 

758 

14 488 

43 667 

466.82 

26.21 

2 051 

8 274 

259.95 



.145 

0 

437 

3 302 

6 925 

46.68 

2.62 

20. 51 

82.74 

26.00 


a Oxide is Cr 2 0 3 ; k y values represent static and moving air (10 times static air values). 



TABLE V. - DIMENSIONLESS TIMES AT CROSSOVER 


AND AT MAXIMUM WEIGHT CHANGE 


Stoichiometric 

constant, 

a 

Maximum- 

specific 

sample 

weight 

change, 

N = a 

Time at 
maximum 
sample 
weight 
change, 

T = /3 

Crossover 

time, 

T o 

T 0 /T 

1.0 

0.3069 

0. 1931 

0. 7969 

4. 126 

1.5 

. 1558 

.1108 

.4506 

4.066 

2.0 

. 09453 

.07213 

.2915 

4.041 

2.5 

.06353 

. 05076 

.2044 

4.028 

o 

CO 

. 04565 

.03768 

. 1515 

4.021 

3.5 

. 03440 

. 02909 

. 1168 

4.016 

4.0 

. 02686 

.02314 

. 09288 

4.013 


TABLE VI. - SUMMARY OF LIMITING VALUES AND 
TRENDS OF THE PARALINE AR EQUATIONS 


Variable 

Limiting value at- 

Time 

zero 

Time 
t = t 

Time 

t=t 0 

Time 
t >> 1 

t 

C 


t 

t 0 

>>1 

AW 



AW 

0 

<0 

W m 



w^ 

>0 

>0 

W 

r 



W 

r 

>0 

(b/a)(A/B) 

T 



P 

T 0 < a- 1 

>>1 

N 



a 

0 

(1/a) - T 

M 



P + (1/b) 

bT 0 

1 + T 

R 



1/b 

fT 0 

1 

dN/dT 

OC 

0 

(T' 1 - b)/a 

-1 




TABLE VII. - DEFINITION OF PARALINE AR PARAMETERS 


AND THEIR CORRESPONDENCE WITH 
TRANSPORT COEFFICIENTS 


Variable 

Definition 

Correspondence 
with transport 
coefficients 

a 

a”* - In (b/a) 


P 

In (b/a) - b' 1 


A 

(/3/F)(AW/a) 2 

a 2 k p /2 

B 

(/Vt)(AW/ a) 

ak v / b 

B/A 

a / AW 

2k y /abk p 

b 2 /a 

P/t 

2k 2 /b 2 k p 

W V 

r, lim 

(b/a)(A/B) 

b 2 kp/2k y 


TABLE VIII. - METHODS FOR CALCULATING 
SYSTEM PARAMETERS A AND B 


Input 

data 

type 

Data given 

Algorithm citation 

1 

AW; ~t 

Table VH 

2 

V k v 

Table Vn 

3 

k p ; AW/t 

Equation (A52) 


k y ; AW/t 

Equation (A57) 

4 

AWj.tj; AW 2 ,t 2 

Equations (A68) and (A69) 


TABLE EX. - SUMMARY OF WEIGHT VALUES AT SPECIAL TIMES 


Weight 

Maximum time, 
t = t and T = 0 

Crossover time, 
t = tQ and T = Tq 

Limiting time, 
t >> 1 and T >> 1 

M 

In (b/a) 

bT 0 

1 + T 

m 

- In (b/a) 
b 

bAtj/B 

A(1 + T)/B 

R 

a’ 2 

aT 0 

1 

W r 

bA/Ba 2 

bAT Q /B 

bA/ab 





Weight change, mg/cm 



Figure 1. - Comparison of parabolic and paralinear kinetics for hypothetical 
alloys with the same parabolic scaling constant (k p = 0. 10 (mg/cm 2 ) 2 /hr) 
but different vaporization constants. Stoichiometric constants, a = 2.0 
and b = 3.0. 
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Paralinear input - Cr-0. 2Y (ref. 4) 


|C 


ic 


c 

c 


c 


/(BLOCK data _ _ 

Dimension hours (ioo) r ppf 1 o i .xxfcg 10 ) . tttmoCo) 
-DIMENSION WW( 100) ,TT (1QO) ,XKP( LOO) , 7QZY ( 1 00 ) 


. COMMON /BLK 1 /HOUR S , QP , XXX , WW, TT 

• 'COMMON /BLK 2 /A SMALL , B SMALL , DENOX , DENMET 

| COMMON/BLK 3 /ITYPEj.LTYPE ,LL , LLX^NUMSET , NUMKV, NUMKP . _ „ 

; OOMMON/CASE 4/XKP,XKV t 

i IlNpUT FOLLOWS , ^ . 

! 1DATA ASMALL ,B SMALL /2 . 16 67 , 3,. 1 6'6 7 / 

• (DATA ITYPE ,LL,LLX/4 ,.l 7, 3/ . , . . . _ „ 

DATA NUMSET./2/ , . ... _ . r . 

DATA QP,/. 5, ,9Q^ . 9.5’/ . L. . 

DATA XXX/71 ,‘9*J5.9., 5 , 719.. 0,/ _ . _ . 

DATA HOURS/., 3 1, ,.6 9, . 8 6.., 1. 1.1 4 ,3 . 06,5 ..83*2Q, 83^ 2.3, 61, 30, 5 6 „ 

*. .5 5. 5 6, .86, 1 1, 116 . 67, 194*4.4, 22,2. 2.2* 2_3 6, 11/ _ . ... 

(DATA WW/2. 50, 3 . 29 , 3. 24, 3. 1.3/ . ______ 

DATA TT/3p,_56, 86.1 L,_2 21,22*236.11/ 

, !end . _ . ... . 


|f.ALS^_j__. _ 

CHROMIUM OXIDE FORMERS . . 

CR.2 0l3l TO. CR03 

_ | |CR- . 2Y TEDMON (REF, 4) BEFORE 


7.19 

5.21 . 1 

& AFTER WMAX I 


_L2_00| JC_. 


1 ATM- 


Descriptors to title, 
output listing 


Input data _ 

. X . . _ 

. 2 _ 

.3 

4 

_ 5_ _ 

6-1 

6-2 

_ . 8 


1A 
2 A 
SA 

4 A 

5 A 


Figure 3. - Typical input data for COREST. Type 4 data; various sets of AWi.tj and AVI^t^ 
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Specific sample weight change, AW, mg/cm 1 


+3.5 


O 
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AWo,t. 



CM 
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+2.5 


+ 2.0 


AWj.tj 

-i 

O 

an 
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□ 

o 


o 


□ 


o 

□ o 
□ 


+1.5 


to 


O Observed 
□ Calculated 


<5 

+ . 5 f^ 


i J l I i 

100 150 200 250 300 

Time, t, hr 

Figure 5. - Specific sample weight change AW observed and cal- 
culated by paralinear analysis (COREST) for Cr 2 03 -forming alloy 

Cr-0. 2Y (ref. 4) at 1200° C. AWj » 2. 50 mg/cm 2 ; AW 2 * 3. 29 
mg/cm 2 ; tj = 30.56 hr ; t 2 ■ 86. 11 hr ; error estimate, 0. 156 
mg/cm 2 . 



69 



Specific sample weight change, AW, mg/cm 



Figure 6. - Specific sample weight change AW observed and calculated by paralinear analysis 
(COREST) for Cr 2 03 -forming alloy Ni-40Cr at 1200° C. Maximum specific sample weight change 
in paralinear oxidation, AW, 3. 24 mg/cm 2 ; time to reach AW, t, 250 hr ; error estimate, 0.243 
mg/cm 2 ; type 1 input. 
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I 



Time, t, hr 


Figure 7. - Specific sample weight change AW observed and calculated by para linear analysis 
(COREST) for Cr^-forming alloy Tophet 30 (2-37) at 1200° C. AWj = 0. 875 mg/cm 2 ; AW 2 = 
-1.061 mg/cm 2 tj = 1200 hr; t 2 = 2050 hr ; error estimate, 0.352 mg/cm 2 



Time, t, hr 

Figure 8. - Specific sample weight change AW observed and calculated by paralinear analysis (COREST) for C^Oj-forming 
alloy Tophet 30 (2-35) at 1150° C. AWi * 1.737 mg/cm 2 ; AWo * 6.330 mg/cm 2 ; ti = 700 hr; to = 2300 hr ; error estimate, 
0.460 mg/cm 2 1 c 
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Sample specific weight change, AW, mg/cm 


IB II II ■■■III ■ 



0 200 400 600 800 1000 1200 1400 1600 1800 2000 

Time, t, hr 

Figure 9. - Specific sample weight change AW observed and calculated by paralinear analy- 
sis (COREST) for Cr^j-forming alloy Tophet 30 (2-42) at 1150° C. AWj = 2.017 mg/cm 2 ; 

A W 2 * 1. 259 mg/cm 2 ; tj = 990 hr; t 2 ■ 1850 hr ; error estimate, 0. 360 mg/cm 2 . 



Time, t, hr 


Figure 10. - CORES T-ca leu lated specific weight of metal consumed W^ for 
predominantly CroOj-forming alloys above 1100° C, assuming paralinear 
oxidation. Stoichiometric constant, a, Z 1667. 
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W m estimates for 1000 hours: 

(1) Assume k p 1 valid only; 

w m, 1000 ' 2. 1667 x 0. 302 x (1000) 112 =■ 20. 7 mg/cm 2 

(2) Assume two-stage parabolic; 

W m, 1000 = 2 - 1667 [°- 302 x < 120)U2 + 0. 129(1000 - 120) 1/2 ] 

= 16. 50 mg/cm 2 

(3) If para linear, k p = 0. 1797 (mg/cm 2 ) 2 /hr ; 

k v = 0. 0140 mg/cm 2 /hr ; W mj 1000 = 32. 61 mg/cm 2 



Figure 11. - Comparison of corrosion estimates (specific weight of metal consumed W m at 1000 hr) if C^C^ -forming Alloy A (ref. 6) were 
analyzed from 300-hour data, assuming parabolic rates. W m = ak^ 2 t 1/2 ; a = 2. 1667 for C^Oj. 



Figure 12. - Comparison of COREST-calculated paralinear and parabolic oxidation results for Cr 2 C ) 3 / 
chromite-forming Alloy A at 1200P C, for same parabolic scaling constant k p . Based on input data from 
reference 6. 
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Parabolic scaling constant, k n , (mg/cm 2 ) 2 /hr 



1200 1150 1100 1050 1000 

Temperature, °C 

Figure 13. - Parabolic scaling constants calculated by COREST, using paralinear 
analysis of literature and this study data, as function of reciprocal of absolute 
temperature. 


Paralinear vaporization rate constant, ky, mg/cm 2 /hr 
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ih 


10 


ri 


10 


-2 


10 


-3 


Curve cr^-derived paralinear vaporization rate 

constants, ky, for - 

A Cr^/chromite alloys (this study and ref. 10) 

B Static air - bulk oxide (ref. 13) 1 ±1 Standard 

C Free vaporization - bulk oxide (ref. 13) j deviation, a 

D Static air - Cro0 3 stripped from oxidized Ni-50Cr 
(ref. 8) 

• Static air - sintered Cr^ (this study) 

Q Alloy A (ref. 6) 

O Cr-0. 2Y (ref. 4) 

A Ni-40Cr (this study) 

D Tophet 30 (this study) 

I Sintered Cr^; 0.2-atm oxygen; airflow, 5.1 to 

29.8 cm/sec (ref. 12) 
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,-4 h- 
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Reciprocal temperature, 1/T, 1/K 


1200 1150 1100 1050 1000 
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Figure 14. - COREST-calculated vaporization rate constants for Cr 2 0 3 as function 
of reciprocal of absolute temperature, assuming paralinear oxidation of Cr^- 
forming alloys. 
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Parabolic scaling constant, k D , <mg/cm 2 ) 2 /hr 
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